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0.  INTRODUCTION 


This  document  summarizes  the  results  of  a  two-year  research  effort 
on  a  number  of  mathematical  and  numerical  issues  related  to  finite  elasto- 
plasticity.  The  Principal  Investigator  of  the  project  was  Professor  J.T. 
Oden  of  The  University  of  Texas. 

The  study  set  out  to  develop  a  general  theory  of  finite  elastoplasticity 
which  makes  use  of  standard  continuum  thermodynamic  arguments,  but  which 
does  not  necessarily  involve  the  assumption  of  existence  of  the  yield 
function,  which  does  not  require  conventional  assumptions  of  convexity 
or  differentiability  of  various  functionals,  which  addresses  and  resolves 
the  issue  of  proper  decomposition  of  elastic  and  plastic  deformation  measures, 
and  which  under  appropriate  additional  assumptions  reduces  to  many  theories 
known  to  be  capable  of  describing  infinitesimal  deformations  of  elasto- 
plastic  solids.  In  addition,  general  forms  of  representation  of  elasto- 
plastic  constitutive  equations  were  sought  which  would  provide  some  measure 
of  the  effects  of  micromechanical  changes  in  the  constitution  of  the  ma¬ 
terial.  Finally,  this  general  theory  was  used  as  the  basis  for  the  con¬ 
struction  of  new  finite  element  methods  for  the  calculation  of  large  elasto- 
plastic  deformations,  together  with  various  algorithms  and  codes  in  order 
to  produce  numerical  solutions  of  representative  problems  in  this  area. 

Positive  results  on  each  one  of  the  objectives  were  obtained.  In 
particular,  the  following  were  accomplished: 

1.  Our  theory  is  a  theory  of  plasticity  that  does  not  necessarily 
involve  the  concept  of  stress:  indeed,  materials  are  fully  characterized 
in  a  way  consistent  with  thermodynamics  by  two  stress  potentials:  the 
free  energy  functional  and  the  generalized  flow  potential. 

2.  The  theory  can  produce  as  special  cases  all  of  the  classical  theories 
of  plasticity  and  hyper-elasticity.  Yet  it  does  not  require  the  introduc¬ 
tion  of  a  yield  function,  nor  does  it  necessarily  involve  assumptions 
of  a  normality  rule  of  convexity. 
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3.  The  notion  of  the  Clark/Rockaf el lar  generalized  sub-differential 
is  used  to  provide  a  broad  generalization  of  classical  potential  theory 
and  make  the  potential  arguments  applicable  to  non-differential  non-convex 
potentials.  Isotropic  function  theory  is  used  to  further  specialize  accept¬ 
able  forms  of  these  functionals. 

4.  A  correct  kinematics  of  finite  elastoplastic  deformations  has 
been  produced,  we  believe  for  the  first  time.  This  involves  the  use  of 
the  polar  decomposition  theorem  and  a  mathematically  correct  definition 
of  the  so-called  stress-relaxed  referenced  state. 

5.  Since  a  correct  kinematics  is  developed  and  potentials  are  used, 
the  formulation  make  unnecessary  the  traditional  difficulties  of  defining 
a  correct,  objective  stress  rate. 

6.  A  new  collection  of  numerical  algorithms  has  been  developed  based 
on  this  theory.  These  have  been  developed  in  connection  with  finite  element 
models  of  finite  elastoplastic  behavior.  A  working  computer  code  has  been 
developed  and  several  example  problems  have  been  solved.  While  the  code 
is  primarily  a  plane  strain  code  which  employs  bilinear  quadrilateral 
elements,  it  has  been  used  to  solve  a  number  of  problems  in  finite  deforma¬ 
tion  plasticity  and  elasticity  and  results  are  consistent  with  those  ob¬ 
tained  by  other  methods  for  special  cases. 

0.2  Publications 

The  following  papers  were  published  as  a  result  of  the  reported  work: 

1.  "Generalized  Potentials  in  Finite  Elastoplasticity ,"  by  S.J.  Kim  and 

J.T.  Oden,  Internat ional  Journal  of  Engineering  Science,  Vol .  22,  No. 

11/12,  1984,  pp.  1235-1257. 

2.  "Generalized  Flow  Potentials  in  Finite  Elastoplasticity  —  II.  Examples," 
by  S.J.  Kim  and  J.T.  Oden,  International  Journal  of  Engineering  Science, 
Vol.  23,  No.  5,  1985,  pp.  515-530. 

3.  "Finite  Element  Analysis  of  a  Class  of  Problems  in  Finite  Elastoplasti¬ 
city  Based  on  the  Thermodynamical  Theory  of  Materials  of  Type-N,"  by  S.J. 
Kim  and  J.T.  Oden,  Computer  Methods  in  Applied  Mechanics  and  Engineering, 
(in  press). 

4.  "KABOD-A  Finite  Element  Program  for  Large  Elastoplastic  Deformations 
Based  on  Generalized  Flow  Potentials,"  by  S.J.  Kim,  TICOM  Report  85-7, 
Austin,  Texas,  1985. 

5.  "Theory  of  Finite  Elastoplasticity,"  by  S.J.  Kim,  TICOM  Report  85-6, 
Austin,  Texas,  1985. 

0.3  Personnel 


The  following  persons  worked  on  this  project: 

1.  Professor  J.  Tinsley  Oden,  Principal  Investigator 
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2.  Seung  Jo  Kim,  Graduate  Research  Assistant,  who  received  his  Ph.D. 
degree  in  May,  1985. 

3,  Tsung-L.  Lin,  Philippe  Devloo,  L.  Demkowicz,  J.  Bass,  and  Luis 
Faria  provided  a  small  portion  of  their  time  during  various  phases  of 
the  project. 

0.4  Technical  Discussion 

A  summary  of  most  of  the  major  results  is  given  in  the  remainder  of 
this  report.  These  results  include  the  development  of  kinematical  relation¬ 
ships,  the  notion  of  generalized  potentials  and  their  thermodynamic  conse¬ 
quences,  representation  results  for  flow  potentials,  finite  element  approxi¬ 
mations,  algorithms,  and  numerical  results,  some  results  on  actual  labora¬ 
tory  experiments  and  their  correlation  with  the  numerical  solutions,  and 
suggestions  for  future  research. 
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CHAPTER  I 

INTRODUCTION 

1.1  Introductory  Comments 

Classical  plasticity  theory  is,  in  many  ways,  a  product 
of  inductive  thought,  growing  as  it  did  from  attempts  to  model 
observed  behavior  of  metals  and  soils  under  loading  histories 
sufficient  to  create  permanent  deformaton.  This  approach  toward 
the  development  of  a  mechanical  theory  is  quite  different  from 
the  deductive  methods  of  modern  continuum  mechanics,  where  the 
framework  for  theories  of  material  behavior  is  derived  in  a  semi- 
axiomatic  way  from  a  small  collection  of  universal  postulates 
(laws  of  physics). 

Over  the  last  fifteen  years,  there  have  been  numerous  attempts 
to  provide  a  thermodynamic  basis  for  a  plasticity  theory  sufficiently 
general  to  accommodate  finite  deformations  but,  at  the  same  time, 
not  inconsistent  with  either  continuum  thermodynamics  or  classical 
plasticity.  While  this  volume  of  literature  is  too  large  to 
be  adequately  referenced  here,  we  mention  as  significant  examples, 
the  works  of  GREEN  and  NAGHDI  (1965],  COLEMAN  and  GURTIN  [1967], 
VALANIS  [1971],  ERINGEN  [1962],  LEE  [1969,  1981],  NEMAT-NASSER 
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[1979,  1982],  HALPHEN  and  NGUYEN  [1975],  BHANDARI  and  ODEN  [1973, 
1975],  RICE  [1971],  HAVNER  and  HILL  [1982],  and  HAVNER  [1982]; 
for  more  recent  work  and  a  more  complete  list  of  references  in 
this  general  area,  the  proceedings  edited  by  LEE  and  MALLET  [1982], 
DESAI  and  GALLAGHER  [1983],  NEMAT-NASSER  et .  al  [1984]  and  WILLAM 
[1984]  can  be  consulted. 

There  are  several  fundamental  features  of  classical  plasticity 
that  are  particularly  difficult  to  generalize  in  a  way  compatible 
with  modern  continuum  mechanics:  the  basic  kinematical  descriptions 
such  as  the  decompositions  of  strain  and  strain-rate  measures 
into  distinct  elastic  and  plastic  parts;  an  appropriate  definition 
of  stress  rate;  the  existence  of  a  yield  function  as  distinct 
from  the  usual  collection  of  constitutive  equations  required 
to  characterize  material  behavior;  so-called  normality  conditions; 
etc.  In  addition,  one  hopes  that  a  theory  which  does  capture 
and  generalize  these  features  will  be  capable  of  modeling  aspects 
of  the  behavior  of  real  materials  which  fall  outside  the  realm 
of  the  classical  theory,  that  it  will  be  consistent  with  continuum 
thermodynamics,  and  that  it  will  be  of  a  form  that  lends  itself 
to  numerical  approximation. 

This  report  presents  the  development  of  a  general  theory 
of  elastoplastic  materials,  which  includes  infinitesimal  elastoplas- 
ticity  as  a  special  case,  and  which  addresses  and  resolves  each  of 


the  difficulties  mentioned  above.  In  addition,  particular  special 
cases  of  the  theory  are  considered  which  provide  meaningful  general¬ 
izations  of  classical  elasticity  and  plasticity  to  cases  of  practical 
importance.  In  addition,  a  study  of  finite  element  approximations 
of  the  governing  equations  is  given,  together  with  new  numerical 
algorithms  and  applications  to  representative  problems. 

1.2  Objectives  and  Scope 

The  first  aim  of  this  study  is  to  develop  a  general  theory 
of  finite  elastoplasticity  which  makes  use  of  standard  continuum 
thermodynamic  arguments,  does  not  necessarily  involve  the  assump¬ 
tion  of  the  existence  of  a  yield  function,  does  not  require  assump¬ 
tions  of  convexity  or  differentiability  of  various  functionals, 
addresses  and  resolves  the  issue  of  proper  decomposition  of  various 
"elastic"  and  "plastic"  deformation  measures  in  a  consistent  way, 
and  which,  under  appropriate  additional  assumptions,  reduces  to 
many  theories  known  to  be  capable  of  describing  infinitesimal 
deformations  of  elastoplastic  solids. 

Secondly,  a  general  form  of  representation  of  the  elastoplastic 


constitutive  relations  is  sought  which  provides  for  some  measure 
of  the  effects  of  micromechanical  changes  in  the  constitution  of 
the  material.  This  is  provided  for  by  regarding  the  principal 


constitutive  functional  as  dependent  on  internal  state  variables 
in  addition  to  appropriate  kinematical  and  thermal  measures. 

Thirdly,  we  wish  to  demonstrate  that  the  theory  lends  itself 
to  numerical  approximation.  Toward  this  goal,  finite  element 
approximations  of  the  governing  equations  are  derived  along  with 
a  collection  of  new  algorithms  for  treating  large-deformation 
elastoplasticity  problems.  A  code  implementing  these  algorithms 
was  developed  for  the  Harris  800  supermini-computer,  and  used 
to  test  the  algorithms  on  a  number  of  representative  example 
problems. 

Finally,  we  wish  to  develop  consistent  finite  element  methods 
for  the  approximation  of  the  equations  and  inequalities  governing 
the  general  theory,  together  with  new  algorithms  for  the  numerical 
analysis  of  representative  problems. 

We  note  that  this  list  of  objectives  is  a  broad  one,  encom¬ 
passing  not  only  the  creation  of  a  new  theory  but  also  new  numerical 
methods  and  actual  applicatons  to  representative  problems.  In 
general,  applications  of  such  a  general  theory  to  concrete  problems 
is  impossible  without  some  information  on  the  specific  form  of 
constitutive  equations  and  accompanying  experimental  data  obtained 
from  tests  on  actual  materials.  To  provide  this  information 
for  the  problem  studied  analytically  and  numerically  here,  we 
construct  slight  generalizations  of  some  flow  potential  functionals 


8 


proposed  by  BODNER  and  PARTOM  {1972,  1975].  These  data  were 
then  used  in  the  determinaton  of  necessary  material  constants. 

This  report  is  divided  into  six  chapters. 

In  Chapter  2,  following  this  Introduction^  mathematical 
and  mechanical  preliminaries  are  given.  The  principal  mathematical 
machinery  needed  for  the  general  techniques  we  employ  is  the 
theory  of  non-convex  optimization  and  generalized  subdifferentials 
advanced  by  CLARKE  (1973,  1976,  1977]  and  ROCKAFELLAR  [1979,  1980]. 
PANAGIOTOPOULOUS  [1982,  1985]  observed  that  these  ideas  had  applicaton 
to  certain  plasticity  problems.  We  outline  the  key  mathematical 
concepts  in  Chapter  2  and  provide  a  brief  review  of  relevant 
ideas  from  non-convex  analysis  in  an  Appendix.  Some  essential 
micromechanical  features  of  materials  relevant  to  plasticity 
theories  for  metals  are  also  reviewed  with  some  interpretations 
of  appropriate  internal  state  variables. 

In  Chapter  3,  the  ideas  and  concepts  introduced  in  Chapter 
2  are  used  in  the  formal  study  of  thermodynamic  restrictions 
on  the  form  of  constitutive  functions  characerizing  a  general 
class  of  materials  which  we  call  "Materials  of  Type  N",  since 
they  involve  a  generalized  "normality"  condition.  A  close  examinaton 
of  the  kinematical  variables  and  a  proposal  of  new  alternative 
kinematical  measures  of  deformaton  and  deformaton  rate  are  also 
presented. 


In  Chapter  4,  general  representations  of  free  energy  and 
flow  potentials  are  developed.  It  is  shown  that  classical  theories 
of  elastoplasticity  can  be  recovered  from  our  theory  as  special 
cases. 

As  an  application  of  the  theory,  specific  forms  of  the  free 
energy  and  flow  potential  functionals  are  presented  which  are 
inspired  by  results  of  BODNER  and  PARTOM  [1972,  1975].  To  demonstrate 
the  feasibility  of  our  theory  and  show  how  to  determine  the  material 
constants  for  our  generalized  Bodner-Partom  material,  actual 
physical  tension  tests  were  performed.  A  numerical  scheme  for 
using  such  sample  experiments  to  determine  material  constants 
is  also  presented  and  discussed. 

The  finite  element  approximation  of  the  governing  equations 
is  taken  up  in  Chapter  5,  where  an  incremental,  total  Lagrangian 
algorithm  is  described.  Specific  applications  are  also  considered, 
including  the  numerical  analysis  of  large-strain  uniaxial  stretching 
and  compressing  of  a  specimen  and  the  crushing  (upsetting)  of 
rectangular  billet  in  a  metal  forming  simulation. 

Among  the  new  and  special  features  of  our  results  are  the 
following: 

Technically,  our  continuum  theory  does  not  require 
a  concept  of  stress.  The  mechanical  response  of  this  class  of 
materials  (excluding  heat  flux)  is  completely  characterized  by 
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two  potentials,  the  free  energy  and  the  flow  potential,  and  only 
subgradients  of  these  functionals  appear  in  the  equations  of 
motion. 

The  potentials  which  play  a  predominant  role  in  our 
formulation  are  true  potentials  only  in  a  very  generalized  sense, 
made  precise  in  Chapter  3;  the  potential  functionals  need  not  be 
convex  nor  need  they  be  differentiable. 

It  is  not  necessary  to  assume  the  existence  of  a  yield 
functon  in  order  to  describe  yielding  and  elastoplastic  deformation; 
however,  if  a  yield  function  is  known  for  a  given  class  of  materials, 
corresponding  flow  potentials  can  be  constructed  in  a  straight¬ 
forward  manner. 

The  various  plasticity  theories  derivable  from  our 
general  formulation  are  not  necessarily  "incremental"  nor  "rate-tpe" 
theories  in  the  usual  sense.  This  fact  opens  the  door  to  several 
new  families  of  numerical  methods  for  solving  finite  deformation 
plasticity  problems.  We  develop  one  such  family  of  new  schemes 
here  and  apply  it  to  representative  examples. 
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CHAPTER  II 


PRELIMINARIES 

In  this  chapter,  several  mathematical  notions  prerequisite  to  our 
study  and  certain  mechanical  concepts  on  plastic  deformation  are  pre¬ 
sented. 


2.1  Mathematical  Preliminaries  and  Non-Convex  Analysis 


We  begin  with  a  summary  of  some  recent  results  on  non-convex 
analysis  which  follows  the  ideas  of  CLARKE  [1973]  and  ROCKAFELLAR 
[1980],  Unless  noted  otherwise,  V  denotes  a  topological  vector 
space,  V*  the  topological  dual  of  V  ,  and  a  duality  pairing 

on  V*  x  V. 


2.1.1  CoYvU.ngo.nX.  and  Tangent  ConeA 

Let  K  he  a  nonempty  subset  of  a  topological  vector  space  V  . 
Then  the  contingent  cone  to  K  at  a  point  u  f  K  is  defined  as  the 
set 

C„(u)  *  lim  sup  (K-u)  (2.1) 

K  0  -  0+  9 


Likewise,  the  tangent  cone  to  K  at  u  is  defined  as  the  set 


Tr(u)  =  lim  inf  -g  [K-u'] 
u*  -►  u 

u*  e  k 
0  -*■  0+ 


-11- 


12 


To  interpret  the  notation  used  in  (2.1)  and  (2.2),  we  use  the 
concept  of  a  limit  superior  (inferior) of  a  multifunction  defined  on 
a  topological  vector  space  given  in  the  Appendix.  Let  T  Be  a  set- 


(2.6) 


U  €  K,3  vk-  v  with  uk  +  ek  vk  f  K}  (2.7) 

A  two-dimensional  case  is  illustrated  in  Fig. 2.1.  Suppose  u 
terminates  at  a  cusp  in  a  non-convex  set  K  ,  as  shown.  The  entire 
plane  can  be  represented  as  the  union  of  four  cones  with  vertex 
u:  BOD,  DOC,  COA,  and  AOB,  with  0  the  terminix  of  u  .  Clearly,  for 
any  point  v  inside  the  cone  BOD  [J  DOC  U COA,  it  is  always  possible 
to  find  a  sequence  of  positive  numbers  {6^}  such  that  u  +  6^  Vr  f  K 
of  any  sequence  v^  -*•  v  .  Outside  of  this  cone  (interior  to  AOB),  it 
is  impossible  to  find  {8^}  for  which  u  +  0^  v^  f  K.  Hence, 

CK(u)  *  BOD  V  DOC  \j  AOB 

Similarly,  pick  a  sequence  u^  -*■  u  where  u^  is  a  sequence  of 

vectors  tracing  out  the  arc  EO  on  K  .  The  legitimate  vectors  v 

with  sequences  v^  -*■  V  such  that  u^  +  0^  v^  £  K  as  0^  -*■  0  will 
be  those  in  the  half  space  BOD  U  DOC  .  Similarly,  for  u^  approach¬ 
ing  u  along  FO ,  we  must  choose  v  in  DOC  U  COA  .  All  other  se¬ 
quences  u^  f  K  ,  u^-*  u  yield  acceptable  v  in  either  of  these 

half  spaces.  Thus,  T  (u)  must  represent  the  intersection: 

TR(u)  -  DOC 


If  K  is  convex,  then 
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2.1.2  Normal  Cone. 

For  K  C  V,  K  4  ft  ,  the  normal  cone  to  K  at  u  is  a  subset 
of  the  dual  space  V*  defined  by 

N  Cu)  =  lu*  f  V*|<u*,v>  <  Q  Vv  fTK(a]}  0.81 

In  two  dimensions,  N^CuL  consists  of  the  vectors  through  u 
which  make  obtuse  angles  with  the  vectors  in  Tj_(u)  »  as  s^own  in 
Fig.  2.2. 

2.1.3  ClaAk-Rodzafiellan.  VeAlva£iveJ> 

Let  F  be  any  extended  real-valued  function  defined  on  V  and 
let  F  be  finite  at  a  point  u  £  V  .  Then  various  types  of  subde¬ 
rivatives  of  F  at  u  can  be  defined  as  follows. 

•  Uppen  Subde/uvative,  The  upper  subderivative  of  F  at  u 
in  direction  v  Is  defined  as 

DtF(u;vl  =  LIm  sup  inf  4[F(u'+0v’ )-a]  (2.9) 

Cu’,aHCu,F(ul),vW  " 

9-  0+ 

The  notation  Lim  sup  inf  is  defined  in  the  Appendix,  and  by  the  no¬ 
tation 


(U;cO  +  (u,F(u) ) 

we  signify  the  convergence  of  a  sequence  (ur,ct)  £  epi  F  to  a  point 
on  the  graph  of  F  : 
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(u'.al  +  tu,F(.u)l<=>(u’ ,ot)  -*■  (u,F(u)), 
a  >  F(u’) 


If  F  is  l.s.c.  Clower  semicontinuous)  at  u  ,  then  (2,8)  re¬ 


duces  to 


D  t  F(u;v) 


Lim  sup 


Cu\F(.uf)  •+•  Cu,F(u))  v*-*-  v 


inf  -ilFCu’+evM-FCu’)] 

f-V  ,,  0 


(2.10) 


•  LoweA.  Suhd&itvcutXve..  The  lower  subderivative  of  F  at  u 
in  direction  v  is  defined  By 


D  +  F (u , v)  =  Lim  inf  sup  -£lF(u’ +  Bv' )  -a] 

(ur,a)  +  (u.F(u)l,  v’+v 
6  -*■  0+ 


(2.11) 


(u',a)  t  (u,F(u))<=>Cu' ,a)  ♦  (u,F(u)) 

a  <  FCu’l 


If  F  is  u.s.c.  at  u  , 


D  1  F(u;v)  =  Lim  inf  sup  -±[F(u'  +6v' )  -F(u' )  1(2. 12) 

Cu’.FCu'))  -*  (u , F (u)  )  ,  v'+  v° 

e  -  o+ 


The  derivatives  DfF(u;v),  D+F(u;v)  are  referred  to  here  as  Clarke- 
Rockafellar  derivatives  (or  C/R-derivatives  for  brevity). 

To  help  understand  the  meaning  of  the  C/R-derivatives,  we  con¬ 
sider  the  two  examples  shown  in  Figs. 2.3  and  2.4. 


Example  1.  Lei  and  f  denote  the  real-valued  functions 

fx2  -  4x  +  4  i  <  1 

fl(x)  H  2 

-x  /2  +  2x  -  1/2  ,  x  >1 


and 


f2(x) 


^ x2  -  4x  +  1 
|-x2/2  +  2x  +  3/2 


x  1  1 
x  >  1 


As  we  observe  in  Fig.  2.3,  is  continuous  at  .:  =  1  but  non- 

dif f erentiable  in  classical  sense  and  f^(x)  is  discontinuous  but 
lower  semicont inuous  at  x  =  1  .  Since  these  functions  are  defined 
on  E  ,  the  definition  of  C/R-derivatives  reduces  to 


Dtf (x;yl 


lim  sup  -i[F(x’  +  9y) 

(x\f(x’)l  Cx,f  (x))  W 

e-o+ 


F(x’)] 


and 

Dtf(x;y)  =  lim  inf  i-[F(x’  +  0y)  -  F(x')] 

Cx\f(x'))  -  (x,f(x))  W 

e  -*•  o+ 

Around  x  =  1  ,  we  have  two  subsequences  of  derivatives  with  y=l 
which  converges  to  1  and  -2  for  both  f^  and  f^  and  we  easily 
conclude  that 

Dtf  ^  (1;  1)  =  1  (=  "sup  { 1 ,  -2  }") 


Dtf  (1; l)  =  1 
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I  D+f  (1;1)  =  -2 

D+f 2 (1; 1)  =  -2 

j  We  remark  that  the  convergence  (x*  ,a)  4-  (x,f(x))  in  the  definition 

of  C/R-derivatives  (2.8)  and  the  definition  of  tangent  cone  (2.2) 

lead  to  the  fact  that  the  epigraph  of  the  function  y  *  D+f(x;y)  is 

j  the  tangent  cone  T^^  ^  (x,f(x))  which  is  shown  in  Fig.  2.3a 

Example  2.  A  better  example  of  C/R-derivatives  can  be  construc- 
2 

ted  in  ]R  .  Consider  the  lower  semicontinuous  function  f  =  f(x,y) 
shown  in  Fig. 2.4.  The  numbers  indicated  in  the  figure  are  intended 

I 

to  mean  the  following: 


The 

slope 

of 

line 

GH 

at 

G 

in 

the 

direction 

V  = 

(1,0) 

is 

+0.5 

The 

slope 

of 

line 

EF 

at 

F 

in 

the 

direction 

V  = 

(1,0) 

is 

-0.5 

The 

slope 

of 

line 

CD 

at 

C 

in 

the 

direction 

V  33 

(1,0) 

is 

-0.8 

The 

slope 

of 

line 

AB 

at 

B 

in 

the 

direction 

V  = 

(1.0) 

is 

+0.2 

! 


Let  us  calculate  C/R-derivative  at  the  origin  u  =  (0,0) 
ection  v  *  (1,0).  Recall  that 

D+f ((0,0) ; (1,0) )  =  Lim  sup  inf 


(u',F(u’))  -  ( (0,0)  ,f  (0, 0) )  vW 
0  -*•  0+  ~  - 

t|(f(u’  +  v')-f(u’)} 


with  dir- 


Along  the  direction  of  v  =  (1,0)  ,  we  choose  two  sequences  approach 
ing  u  =  (0,0)  from  either  positive  or  negative  x-axis: 


2.4  The  upper  (lower)  subderivat 
continuous  function  F(u). 
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Lim  inf  Df(u’,v')  taken  from  the  positive  side  of  x-axis  will  be  a 
v'  v 

sequence  of  slopes  along  EF  and  Lim  inf  Df(u',v*)  taken  from  the 

v'  -►  v 

negative  side  of  x-axis  will  be  the  slopes  along  DC,  where  Df(u',v') 

*  i[f(u'  +  8v’)  -  f(u')J.  Next  we  take  lim  sup 

Cu'.f(u'))  4-  ((0,0),f(0,0)) 

of  the  sequences  of  slopes  and  then  finally  the  slope  at  F  along  the 
curve  EF  as  Dtf ((0,0) ; (1,0)) .  Similarly,  we  can  get 
D4f ((0,0) ; (1,0))  *  the  slope  at  B  along  AB. 

Quantitatively 

D+f ((0,0) ; (1,0))  -  -.5 

and 

D4-f  ((0,0)  ;  (1,0))  *  +.2 

2.1. A  Generalized  SubcU^exziitiaZi 

The  subdifferential  3F(u)  of  a  function  F  at  u  is  a  well  known 
concept  in  convex  analysis  (see  e.g.,  EKERLAND  and  TEMAM  [1976]  or 
ODEN  [1985]).  By  using  the  subderivatives  defined  in  (2.8)  and  (2.10), 
we  define  the  generalized  subdifferential  of  F:V  +  E,  at  a  point  u 
where  F(u)  is  finite,  as  the  set 

3F(u)  *  {u*  (  V*  |  <u*,v>  £D+F(u;v)  Vvf  V}  (2.13) 

Now  we  list  two  useful  theorems  due  to  ROCKAFELLAR  [ 1980] . 

Theorem  2.1.  Let  F  be  any  extended  real-valued  function  on  V, 
and  let  u  be  any  point  at  which  F  is  finite.  Then  3F  is  a  weak*- 
closed  convex  subset  of  V*  and 


3F(u)  -  (u*  (  V*  j  (u*,-l)  (  NepiF(u*F(u))}  (2* 3.4) 

If  DtF(u;0)  =  then  3F(u)  is  empty,  but  otherwise  ^F(u)  is 
nonempty  and 

D+F(u;v)  =»  sup  {<u*,v>  |  u*  f  lF(u)}for  all  v  f  V  (2.15) 


Theorem  2.2.  If  F  is  a  convex  function  on  V,  then  3F(u)  agrees 
with  the  subgradient  set  in  the  sense  of  convex  analysis: 
lF(u)  “  3F(u) 

“  lu*  f  V*  |  <u*,v)  £  F'(u;v),Vv  6  V} 


(u*f  V*  |  ^v— u,v*}><_  F(u)-F(u)  V  v  (  V} 


(2.16) 


Here  F'  (u;v)  ■  lim  (F(u+tv)  -  F(u))/t  is  called  the  one-sided  direc- 

t-0+ 


tional  derivative  which  exists  for  all  v  when  F  is  convex  (although  it 
may  be  infinite). 

Remark.  If  F(u)  is  a  characteristic  function  with  respect  to 


a  set  K,  i.e.,  if 


F(u)  -  ipK(u) 


0  if  u(  K 
+«  if  u^K 


then 


3F(u)  =  {u*  £  V*  |  <  u*,v )  <_  DtH^ujv)  y  v  f  V)  (2.17) 

-  tyu) 


This  fact  can  be  more  easily  visualized  in  the  case  of  convex  F,  i.e.. 


3^k(u)  =  {u*  £  V*  (<u*,  v-u  > 

— i^K ( v >  “  V  v  f  K} 

■  {u*  £  V*  |  (  u*,  v-u  y  <  0, V  v  £  K; 


Vu) 


since  i>R(v)  =  ^R(u)  =  0. 


2.2  Mechanical  Preliminaries 


2.2.1  Classical  Continuum  Plasticity 

Classical  plasticity  theory  rests  on  the  assumption  of  the 

3  3 

existence  of  a  convex  function  F  :  1M  *  ]R  x  ®  -*■  [0,-®  )  of  the 

stress  tensor  a,  called  the  yield  function  of  the  material,  which 

has  the  property  that  plastic  flow  at  a  particle  X  of  the  material  is 

signaled  whenever  F(a  (X))  =  0  and  ~  :ct>0;  otherwise  the  deformation 

■>»  dO  —  — 

T 

at  X  is  elastic  (where  A:B  =  tr  A  B): 


F(  a  (X))  <  0 


F(  a  (X)  )*  0  and  —  :  a  <0 

~  ~  da  ~  J 


elastic  deformation 


(2.18) 


F(  o  (X))*  0  and  —  :  a  >  0  =>>  plastic  flow 

~  ~  da  — 


The  only  stress  states  admissible  in  such  theories  are  those  for 
which  F(o)  <0  or,  equivalently,  those  stresses  which  belong  to  the 


convex  set 


K  «{a£  3M  !  F(  o)  <0}  C2.19) 

In  general,  the  yield  function  may  itain  several  parameters 
characterizing  the  state  of  material,  but  during  plastic  deformation 


these  parameters  are  understood  to  change  in  such  a  way  that  F 
remains  equal  to  zero  while  plastic  flow  continues.  These  material 
parameters  will  also  be  expected  to  vary  with  temperature  or  other 
thermodynamic  variables.  The  general  yield  condition  relates  only 
to  the  state  of  stress  at  a  material  particle,  and,  irrespective 
of  whether  the  mechanical  response  is  elastic  or  plastic,  does  not 
depend  on  the  stress  gradients. 

It  is  further  assumed  in  the  classical  theory  of  plasticity 
that  yield  is  unaffected  by  hydrostatic  stress,  a  situation  generally 
in  accord  with  experiment  for  moderate  stress  levels.  With  this 
assumption  it  is  possible  to  represent  the  yield  function  for  an 
isotropic  material  as  a  function  of  the  invariants  of  the  stress 
deviator  in  the  form 

9(a)  =■  F  (J2,J3)  (2.20) 

Where  the  invariants 

J2  *  1/2  tr  (o')2,  *  det  a'  with  a'*  a- 1/3  tr  a  (2.21) 

(More  detailed  discussions  of  representations  of  response  functions 
are  given  in  Chapter  IV.) 

The  infinitesimal  strain  tensor  e  is  representable  as  the 

e  P 

sum  of  an  elastic  part  e  and  a  plastic  strain  e  ,  and  its  time 


rate-of-change  is 


It  is  meaningful  to  assume  the  existence  of  a  plastic  flow 
potential  ^  :  1M  -*■  R  which  is  convex  and  l.s.c.  and  which  has  the 


property  that 


i P(o) 


(2.22 


In  particular,  the  indicator  function  \pv  of  the  set  K  may  define  a 


specific  flow  potential  as  follows: 


V°J 


f  0  If  Of  K 

V  +  10  i f  0 


(2.2  L 


Note  that  ip  is  l.s.c.  on  IM  and  that  ip  is  convex  if  F  is  convex. 
K  K 

From  (2.23)  and  the  definition  of  the  subdifferential 


t  n*K  (a) 


(2.2! 


for  some  particular  stress  a  ,  and  this  implies  that 

<^e  P,  a-  a>  >0  \f  af  K  (2.2( 

This  result,  of  course,  is  the  classical  normality  condition  which 
extablishes  that  the  strain  rate  is  normal  to  the  yield  surface  or 
lies  in  the  normal  cone  of  the  yield  surface  at  corners.  (See  Fig. 
2.3  ). 


2.2.2  mcA.omzchaiu.ccU  Plcu>£icU.y 

It  is  now  a  well-known  fact  that  the  permanent  deformation  of 
crystalline  material  is  caused  by  microscopic  defects  in  the 
crystalline  structures.  A  brief  review  of  some  features  of  these 
underlying  physical  phenomena  is  useful  at  this  point.  For  a 
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detailed  treatment,  COTTRELL  [ 1958 ] ,  WEERTMAN  and  WEERTMAN  [1964], 
and  LARDNER  [1974]  can  be  consulted. 

When  most  materials  solidify  from  a  melt,  the  constituent  atoms 
arrange  themselves  in  a  highly  regular  periodic  array  called  a 
crystal.  They  do  this  simply  because  the  total  energy,  kinetic  and 
|  potential,  of  the  individual  atoms,  is  less  when  the  atoms  are 

arranged  in  a  coordinated  manner  than  when  they  are  randomly 
distributed.  This  energy  minimum  changes  for  different  atomic 
constituents,  so  that  the  stable  crystal  structure  formed  during 
solidification  vary  for  each  material.  The  exact  nature  of  the 
structure  is  determined  by  the  strength  and  directionality  of  the 
interatomic  forces. 

A  high  degree  of  geometrical  perfection  in  crystals  is  rarely 
found  in  either  natural  or  synthetic  crystals.  In  all  real 
materials,  crystals  contain  a  small  number  of  atomic  irregular¬ 
ities  or  imperfections  which  do  little  to  change  the  normal  atomic 
arrangement  but  play  an  important  role  in  modifying  its  physical 
properties.  The  defects  of  most  importance  in  determining  physical 
behavior  are  dislocations,  point  defects,  and  the  structure  of  grain 
and  phase  boundaries. 

The  dislocations  act  as  the  carriers  of  plastic  deformation 
through  the  motion  of  slip,  and  produce  most  of  the  noticeable  changes 
of  the  internal  mechanical  state  during  the  deformation.  Point 


defects,  such  as  vacancies  and  solute  atoms  of  a  different  species 
usually  introduce  a  strong  viscous  effect  into  the  deformation 
behavior.  Finally,  the  structure  of  the  grain  and  phase  boundaries 
not  only  can  increase  the  resistance  against  plastic  flow  but  they 
can  also  be  additional  sources  of  viscous  behavior. 

For  a  perfect  single  crystalline  material,  one  would  expect 
elastic-perfectly-plastic  response  because  when  the  shearing  stress 
reaches  some  critical  value  it  causes  a  continuous  slipping  motion  as 
there  are  no  obstacles  which  may  interrupt  or  stop  the  shearing  motion 
But,  as  mentioned  earlier,  most  materials  do  have  barriers  like 
grain  boundaries,  Frank  networks,  stationary  substructures,  etc.  To 
break  these  barriers,  the  critical  stress  (or,  phenomenologically,  the 
yield  stress)  must  be  attained.  This  increased  stress  produces  more 
barriers  and  again,  an  increased  stress  is  required  to  break  dis¬ 
location  barriers,  and  so  on.  This  phenomenon  is  called  work¬ 
hardening. 

Generally,  in  single  crystalline  models,  we  have  four  distinct 
stages  in  the  stress-strain  curve, as  seen  in  Fig.  2.6.  After  an 
elastic  stage,  one  often  observes  an  "easy-glide"  stage  (Stage  I  in 
the  figure)  during  which  all  the  free  (mobile)  dislocations  move  and 
a  large  amount  of  plastic  strain  is  realized.  Stage  II  is  called  the 
work-hardening  stage,  and  unit  dislocations  are  generated  from,  say, 
Frank-Read  sources  and  iuteractions  between  dislocations  making 
barriers  like  Lomer-Cottrell  locks.  In  Stage  III,  breakdowns  of 


these  barriers  may  occur,  giving  very  low  hardening  with  some  plastic 
strain.  For  more  detailed  descriptions  of  these  concepts  ,  see 
DIETER  [1976]  or  WILKOV  [1983]. 

In  polycrystalline  materials,  the  easy  glide  stage  is  not 
observed  as  frequently  as  in  the  single  crystal  since  the  movements  of 
dislocations  are  stopped  more  readily  by  a  large  density  of  built-in 
obstacles  due  to  the  nature  of  the  grain  boundaries.  Therefore  the 
usual  stress-strain  curve  (Fig.  2.7),  i.e.,  a  linear  elastic  region 
and  a  strain  hardening  region,  is  obtained.  In  fact,  the  easy  glide 
stage  is  included  in  the  elastic  region  since  the  flow  of  dislocation 
and  plastic  strains  are  ignorable. 

2.2.3  Internal  State.  \)aru.a.bleJ> 

From  the  previous  section,  we  recognize  that  changes  of 
internal  structures  in  materials  causes  permanent  deformations.  A 
key  question  is  how  can  one  utilize  this  micromechanical  information 
in  formulating  a  continuum  theory  of  elasto-plasticity .  One  approach, 
which  has  provided  for  some  progress  toward  a  general  theory,  involves 
the  use  of  continuum  theories  which  employ  "hidden  variables"  or 
"internal  state  variables"  as  a  measure  of  physical  changes  in  the 
microstructure.  Some  of  these  changes  can  be  those  ordinarily 
associated  with  plastic  deformation.  These  so-called  internal  state 
variables  can  provide  measures  of  slip  rearrangements  of  crystall¬ 
ographic  planes  through  dislocation  motions,  inelastic  behavior 
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arising  from  twinning  in  crystals,  grain  boundary  sliding,  and 

stress-induced  phase  transformations. 

The  idea  of  developing  a  continuum  theory  of  thermodynamics 

using  internal  state  variables  was  actually  proposed  by  MEIXNER 

[1953]  and  independently  by  BIOT  [1954].  A  general  thermomechanical 
theory  of  materials  with  internal  state  variables  was  developed  by 
COLEMAN  and  GURTIN  [1967]  and  used  in  attempts  by  PERZYNA  [1963], 
PERZYNA  and  WOJNO  [1968],  ODEN  and  BHANDARI  [1973],  KRONER  [1962], 
KRONER  and  TEODOSIU  [1974],  RICE  J 1971,  1975]  and  others  to  develop 
a  plasticity  theory.  Among  the  advantages  of  such  an  approach  are 
that  it  simplifies  the  modeling  of  history  effects  in  elastic-plastic 
materials  since  they  are,  in  essence,  accounted  for  by  an  evolution 
equation  for  the  internal  state  variables.  There  are  three  categories 
of  classifications  of  internal  state  variables. 

The  first  employs  the  classical  theory  of  elastop last ic ity  and 
its  variations,  e.g.,  MROZ  [1973,  1981]  and  NECAS  and  HLAVACEK  [1981]. 

In  this  approach,  the  hardness  variables  in  the  flow  theory  and  the 
back  stress  in  the  kinematic  hardening  are  the  internal  state  variables. 

In  the  second  approach, a  direct  application  of  micromechanical 
structure  is  modeled,  e.g.,  KRONER  [1962],  KRONER  and  TEODOSIU  [1974], 
RICE  [1971],  ACHENBACH,  MULLER  and  WILMANSKI  [1981]  and  WENG  [1980, 

1981]  Various  measures  of  dislocation  density,  degree  of  impurity, 
and  twinning  modes  can  be  identified  as  internal  state  variables. 


The  third  approach  employs  a  combination  of  the  ideas  of  internal 
state  variables  and  the  phenomenological  behavior  as  derived  from 
experiments,  e.g.,  HART  [1976,  1979],  BODNER  and  PARTOM  [1972, 

1975],  etc.  This  approach  combines  features  of  the  other  two.  More 
complex  constitutive  equations  with  more  variables  which  are  drawn 

from  the  phenomenological  results  are  involved.  A  simple  version  in 
this  category  may  be  an  elastoplasticity  theory  (ARAVAS  [1985])  with 
micro-void  implementation  in  which  a  classical  yield  function  (or 
flow  potential),  e.g.,  von  Mises,  is  modified  by  introducing  the 
volume  fraction  of  micro-voids  in  the  material.  This  particular  theory, 
which  may  find  rare  applications  for  engineering  materials  undergoing 
large  deformations  (see  FUNG,  BURNS  and  LIND  [ 1973 ]),  results  in  a 
pressure  dependent  plasticity  theory  which  may  be  capable  of  modeling 
a  class  of  porous  materials.  Later  we  utilize  this  third  approach 
in  an  application  of  the  theory  of  Materials  of  Type  N. 


CHAPTER  III 

THE  MATERIALS  OF  TYPE  N 


In  this  chapter,  a  general  continuum  theory  is  developed  for 
a  hypothetical  class  of  materials  which  undergo  large  plastic 
deformations . 

3.1  Kinematical  Considerations  in  Finite  Deformations 


While  somewhat  difficult  to  define  with  flawless  mathematical 
precision,  the  basic  idea  of  a  finite  "plastic"  deformation 
of  a  material  is  heuristical ly  clear:  it  is  the  "irrecoverable" 
part  of  the  deformation  of  a  material  subjected  to  a  loading 
cycle.  One  imagines  that  a  material  in  a  reference  stress 
state  Qq  at  a  particle  X  is  subjected  to  motions  which  carry 
the  stress  at  this  particle  through  a  history  which  eventually 
returns  the  stress  to  the  original  state  o^.  If,  at  the  conclusion 
of  this  stress  cycle,  the  local  state  of  deformation  at  X,  however 
one  chooses  to  measure  it,  differs  from  what  it  was  before  the 
stress  cycle,  then  a  portion  of  the  deformation  was  not  "recovered" 
and  this  is  dubbed  the  "plastic"  deformation.  It  is  clear  that 
these  ideas  are  local  in  character;  they  may  have  meaning  only 
in  a  local  neighborhood  of  a  material  particle  or  for  bodies 
in  states  of  homogeneous  deformation.  Having  these  ideas  in  our 
mind,  we  offer  in  this  section  alternative  decompositions  of 
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the  deformaton  gradient  to  the  multiplicative  and  the  additive 
ones  proposed  by  LEE  [1969,  1981]  and  NEMAT-NASSER  [1979,  1982], 
respectively. 


Let  us  consider  the  motion  of  material  body  B  relative 
to  a  fixed  configuration  Cq  C  IR^  (N  _<  3),  which  is  defined 
by  the  map  B  ♦  1RN,  X  =  Kq(X),  where  X  is  a  material  particle. 
The  spatial  position  x  of  a  particle  X  at  time  t  is  then  given 
by  a  relation  of  the  type 


x  =  x(X.  t)  (3.1) 

with  X  £  <q(B),  t  _>  0,  and  x  a  continuous  invertible  map  from 
N 

Cq  into  1R  .  The  deformation  gradient  tensor  F  at  X  at  time 
t  is  defined  by 

SX 


F 


(3.2) 


Let  N  (X)  denote  a  small  material  neighborhood  of  particle 
X.  The  motion  of  the  body  carries  N  (X)  from  the  reference 
configuration  Cq  to  the  current  configuration  C^.  Let  the  Cauchy 
stress  o  at  any  particle  A  £  N  (X)  in  Cq  be  denoted  o (A,  0)  with 
o(A,  0)  *  Oq(X)  +  w(AX)  V  A  6  N  (X) 


where  AX  =  A  -  X.  If  o  is  continuous  at  X,  then 


lim 

IMbo 


|"(  fix) 


||  *11  denotes  the  Euclidean  norms.  We  shall  refer  to  °q(x)  as 
the  initial  stress  at  particle  X.  For  _  .plicity,  we  omit  other 
variables  (such  as  temperature,  etc.)  that  could  also  be  listed 
in  defining  the  "initial  state"  of  the  material. 
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During  an  interval  of  time  [0,t],  t  >  0,  the  stress  at 

particles  in  W  (X)  are  part  of  the  stress  history 

Ht(A)  *  {c(A,  t )  |  A  fW  (X),  0  <_  t  <_  t} 

and  the  configurations  of  W(X)  are  denoted  x(  N  (X),  x  ), 

0  1  r  1 

In  addition  to  the  actual  stress  history  H^,  we  consider 
£ 

any  stress  history  H^,  corresponding  to  a  relaxation  of  the 


stress  at  X,  such  that 


H*(A)  C  -t^*  -t  f  H-  0 


where  is  the  family  of  all  stress  histories  terminating 


at  !o: 


H*(A)  e  H  0  ->  0( A,  t)  =  o0(A,  0) 


For  p  the  values  of  a  continuous  map  of  N  (X)  into  ]R  , 

we  denote  by  any  configuration  of  B  for  which  the  stress 

history  at  A  f  N  (X)  is  Ht(A).  Thus,  the  introduction  of  a 

(possibly  unattained)  configuration  provides  for  the  familiar 

device  of  comparing  the  geometries  of  the  body  in  with  Cq 

to  define  plastic  deformation.  A  typical  illustration  of  this 

is  given  in  Fig.  3.1,  wherein  N  (X)  is  shown  in  Cq,  C^,  and 

C  together  with  differential  vectors  dX,  dx,  and  dp  [cf.  LEE 
P  ... 

[1981,  p.  862]]. 

We  shall  now  proceed  to  construct  decompositions  of  the 


total  deformation  into  elastic  and  plastic  parts. 
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3.1.1  NEMAT-NASSER'4  PecompoAxL^con 

We  first  consider  NEMAT-NASSER' s  additive  decomposition 
which  can  be  obtained  by  manipulation  of  a  displacement  field 
as  follows  [cf.,  NEMAT-NASSER  [1979,  p.  166)]: 

dx  -dX  =  (F  -  I )dX  *  du  (3.3) 

dp  -  dX  -  (FP  -  I)dX  -  duP  (3.4) 

If  we  demand  that  a  single-valued  displacement  field  must  be 
realized  in  reaching  the  current  configuration,  we  are  led  to 
a  definition  of  an  elastic  deformation  gradient  given  by 

( u(X)  +  du)  -  [ u(X)  +  dup]  -  (Fe  -  I)dX  (3.5) 

We  next  use  equations  (3. 3), (3. 4)  and  (3.5)  to  eliminate  du  and 
dup  and  obtain 

F  =  Fe  +  FP  -  I  (3.6) 


In  other  words, 

F*dX  -  (du  -  duP  +  dX)  (3.7) 

and 


FpdX  -  dp  =  (duP  +  dX)  (3.8) 

Since  the  material  derivative  of  x  is  defined  by  x=3x(x,t)/3t=u. 


jl  i i  * -1 

3x  3  t  3~t\3  X/3  ::  =  “ 


(3.9) 


and,  from  equation  (3.6) 


F  -  Fe  ♦  FP 


(3.10) 


We  have 


L  *  FeF_1  ♦  FPF_1  «  L6  +  LP 


(3.11) 


Finally,  the  symmetric  part  of  L  is  given  by  the  sum, 

D  =  i(L  +  LT)  =  i(L6  +  LeT)  +  i(LP  +  LpT)  (3.12) 

=  De  +  DP 

where  D®  =  FeF_1|gymDP  =  (3.13) 


and  I sym  denotes  the  symmetric  part  of  the  tensor. 


3.1.2  LEE 'a  Ve.compoAlti.on 

Lee's  decomposition  starts  "from  the  chain  rule"  [see  LUBARDA 
and  LEE  [1981,  p. 4 ] : 


F  =  &  =  h.  ,  f®Fp 

3X  3  p  3  X  :  5 


Then,  by  using  equations  (3.11)  and  (3.14), 

L  =  FF_1  *  (F^FP  ♦  F'eFP)FP"1Fe"1 
=  feFe_1  +  FeFPFP_1Fe-1 


-e  +  peEp-e-l 


(3.14) 


(3.15) 


We  recognize  from  (3.5)  that  the  rate-of-strain  measure  D  could 
not  be  decomposed  additively,  as  Lee  noted.  In  metals,  elastic 
strain  rates  may  be  small,  and  equation  (3.12)  can  be  approximated 
by 

L  :  L®  +  LP  and  D  =  D®  +  DP  (3.16) 

But  this  additive  decomposition  is,  in  general,  invalid  for  the 
(perhaps  unusual)  case  of  finite  elastic  deformations. 

Although  we  want  a  strain  rate  which  is  independent  of  the 


rotation  effects,  the  previous  decompositions  do  not  render  a  full 


stretch-dependent  strain  rate  since,  e.g.,  Dp  =  FKF  |  m  and  F 
contains  the  rotation  effect  as  F  =  RU  =  VR.  Furthermore,  the  previous 
decompositions  pose  a  question  about  the  conceptual  validity  of 
the  differential  operation  dp/dx.  This  fact  may  generate  an  overdeter¬ 
mined  system  of  equations.  The  determinacy  of  a  system  of  elasto- 
plastic  equations  will  be  discussed  in  a  later  section  of  this 
chapter  after  all  the  field  equations  and  the  constitutive  equations 
are  established  through  the  theory  of  materials  of  type  N. 

To  cure  the  difficulties  mentioned  above,  we  will  search  for 
new  kinematical  variables. 

3.1.3  klZ&inaXA.'JZ  Vzcompo Auction  HeXhocLh 

Consider  again  a  particle  A  in  the  neighborhood  W ( X )  of  X, 

in  the  reference  configuration  (see  Fig.  3.2  or  Fig.  3.3).  The 
position  vector  of  point  A  relative  to  the  origin  of  the  fixed 
spatial  reference  frame  is  denoted 

0Ar  *  X  +  AX 
C0  ' 

According  to  the  polar  decomposition  theorem  (see  e.g.  GURTIN 
[1981]),  the  deformation  gradient  F  can  be  represented  as  the  composi¬ 
tion, 

F  =  RU  =  VR  (3.17) 

where  R  is  a  positive  definite  orthogonal  rotation  tensor,  U  is 
the  right  stretch  tensor  of  F,  and  V  is  the  left  stretch  tensor  of  F. 


A)  Right  Polar  Decomposition 

The  location  of  A  in  the  rotation  free  current  configuration 
C  ,  as  in  Fig.  3.2,  is 

0AC  =  x  +RTAx  *  x(X,t)  +  U(X,t)AX  +  W  (X,  t,  AX) 

t_  (3.13) 

where 

lim  flTTTT  Hwi<x*  t,  AX)  ||  =  0  (3.19) 

Haxim  a  ' 

Now  let  us  consider  an  intermediate  configuration  C^_  which 
may  correspond  to  a  "rotation  free"  state  at  the  initial  stress 
level  Oq.  Then, 

0A_  -  X(X,t)  +  UP(X,t)AX  +  W  (X,  t,  AX)  (3.20) 

P- . 2  ' 

and  W2  has  the  same  asymptotic  behavior  with  respect  to  ||  AX  ||  as 
does  Wj . 

e 

We  next  introduce  a  second  order  tensor  U  which  represents 

the  elastic  stretch  tensor  and  is  defined  by 

Ue(X,t)AX  -  AX  »  0AC  -  0AC 

t-  p- 

=  [U(X,t)  -  UP(X,t)]AX 

♦  W}(X,  t,  AX)  -  W  (X,  t,  AX) 

(3.21 ) 

Thus,  in  the  limit  as  |AX|  *  0,  we  have 

U  =  Ue  ♦  UP  -  1  (3.22) 


and  the  relations  (3.17>  and  (.3.11)  give 


F  =  RU  *=  RU6  +  RUP  -  R 
L  -  F  F_1 


(3.23) 


=  R  RT  +  R(UeU_1  +  UPU_1 )RT  (3.24) 


e  p 

Therefore,  D  and  Dr  can  be  obtained  as 

D®  -  R  Ue  if1  RT|  (3.25) 

.  ^  _  sym 

DP  *  R  UP  U_1  RT|  (3.26) 

B)  Left  Polar  Decomposition 

The  left  pclar  decomposition  cannot  be  constructed  as  easily 
and  simply  as  the  right,  but  a  formal  manipulation  of  the  key  tensors 
is  still  possible.  In  this  case,  the  material  is  first  rotated 
as  in  Fig.  3.3.  Therefore,  the  location  of  A  in  the  rotated  reference 
configuration  Cq+  is 


0AC  =  X  +  R  AX  (3.27) 

whereas  the  locations  of  A  in  the  current  configuration  Ct  and  in 

the  rotated  stress-free  configuration  C  are 

P+ 

OA^  *  x  +  Ax 
t 

=  x(X,t )  +  V(X,t)RAX  +  W3(X,  t:  AX) 

(3.28) 

4 

and 

0A„  =  x(X, t )  +  VP(X,  t  )RAX  +  W.(X,  t  AX)  (3.29)  j 

P  i 


1 


Here  and  have  the  same  asymptotic  property  as  do  and  W 

In  a  way  similar  to  (3.21),  a  second-order  elastic  (left) 
0 

stretch  tensor  V  is  defined  by 

V6(X,  t )RAX  -  RAX  =  0AC  -  0AC 

t  p+ 

=  [V(X,  t)  -  VP(x ,  t ) ]RAX 

+  W  (X,  t,  AX)  -  W  (X,  t,  AX) 

'  '  '  -  '  '  (3. 


And  ||  AX  ||  ♦  0  gives 


V  =  V6  +  VP  -  1  (3. 

or 

F  “  VR  =  V6  R  +  VP  R  -  R  (3. 

Finally,  from  the  definitions  of  L  and  D, 

*  -1 
L  =  F  F 

=  V6  V  ^  +  vP  +  V  R  R^  V-^ 

=  ve  v"1  +  VP  v_1  +  gj(Ve,  R)  ♦  g2(VP,  R) 

(3. 


where,  for  example, 

gj(VeR)  =  VeR  RTV-1  and  g2(VP,  R)  =  VPR  RTV_1  -  R  rV 

(3. 


Then  one  could  define, 
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r 


and 

Dp  =  (Ve  V_1  +  g_(VP,  R) ) |  (3.56) 

It  is  noted  here  that  unlike  NEMAT-NASSER  [1979,  1982]  or 
LEE  [1969,  1981],  neither  Ve  nor  Vp  need  be  a  function  of  the  gradient 
of  an  "elastic"  or  "plastic"  displacement  or  position  vector  field. 
However,  if  u  is  the  particle  displacement  field,  we  always  have, 

F  =  1  +  Vu  (3.37) 

It  has  recently  been  pointed  out  to  us  that  our  decomposition 
through  (3.17)  is  similar  to  that  proposed  by  SIMO  and  MARSDEX  [1984]. 

3.2  Materials  of  Type  N 


The  thermomechanical  behavior  of  the  body  is  governed  by  the 
principles  of  conservations  of  mass  and  energy,  balances  of  the 
linear  and  the  angular  momenta,  and  the  law  of  entropy  production. 
Local  forms  of  these  principles  can  be  written  as  follows: 
•Conservation  of  Mass 

p  det  F  =  pQ  (3.38) 

•Balance  of  Linear  and  Angular  Momenta 

div  o  +  p  b  =  px  (3.39) 

o  =  aT  (3.40) 


•Conservation  of  Energy 

P  (f  =  (a :  L)  -  div  q  +  pr 


(3.41) 


(3.42) 


'Clausius  -  Duhem  Inequality 

•  q  r 

Pn  +  div  j  -  0  _>  0 

Here  P  is  the  mass  density,  PQ  the  mass  density  in  the  reference 

configuration,  a  is  the  Cauchy  Stress,  b  the  body  force  per  unit 

mass,  6  the  specific  internal  energy,  q  the  heat  flux  vector,  r 

the  heat  supply  per  unit  mass  per  unit  time,  n  the  specific 

entropy,  and  9  the  absolute  temperature. 

It  is  customary  to  introduce  the  free  energy  density  function 

instMH  of  the  internal  energy  density  £  by 

♦  *  €  -  n0  (3.43) 

Now  we  are  at  the  place  to  define  a  general  class  of  materials 

which  satisfies  certain  properties. 

Let  us  first  introduce  sets  SCIR3xIR3,  A  C  K3  x  ]R  3 

and  W  «  S  x  A. 

A  material  is  said  to  be  of  Type  N  if  and  only  if  it  is 
characterized  by  constitutive  equations, 

<t>  *  *  (E,0 ,g,a) 
o  =  l  (E,0 ,g,a ) 

.  (3.44) 

n  *  N  (E,0  ,g,a  ) 

q  =  0.  (E,e  ,g,a) 

and  there  exists  a  potential  ip:W  *  (-oo,®).  Which  has  a  nonempty 
generalized  subdifferential  such  that 

(P,-a)<f  3iKo,A),  V(o,M  f  W 


(3.43) 
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Here  g  *  grad  9  and  is  a  tensor  referred  to  as  an  internal  state 
variable.  The  elastic  strain  measure  E  and  the  plastic  strain 
measure  P  were  used  without  specific  relations  with  the  kinematical 
variables.  We  will  investigate  these  for  all  four  cases  which 
were  introduced  in  Section  3.1. 

From  the  result  (2.13),  the  relation  (3.45)  implies  the 
inequality 

<P.o*>s  ♦  <-d,A*>A<_  D+4i((o  ,A)  ;(a*,A*) ) 

V(a*,A*)  £  U  (3.46) 

Where  <  v)s  and  <  • ,  •  >  denote  duality  paring  on  S*  x  S  and 
A*  x  A,  respectively. 

Geometrically,  we  have  the  normality  (hence.  Type  N)  condition, 
(( p »  “<*).  'l)f  Nepi^°’^’  (3.47) 


This  relation  is  illustrated  graphically  in  Fig.  3.4  with 
U  ■  (a, A)  and  F  »  ty. 


Furthermore,  if  the  potential  ij,  is  differentiable  at  (a, A) 


then  N  ..at  (o»A)  has  a  single  element.  So, 
epnj>  ,  _ 


(P,  -a) 


d(n,A) 


or 


and 


p  =  Jjk 


-/ 


(3.48) 


Nepi  fr(U2  ,F(U,)) 


Ncticing  that  none  of  the  constitutive  variables  in  the  assump- 

•  • 

tion  of  materials  of  Type  N  are  functions  of  9  and  g,  we  can  vary 
the  variables  9,  g  without  changing  the  other  variables.  Then  we 
will  obtain  the  following  equalities 


n  -  N(E,0,g,&)  -  - 

0  “  ¥ 

3g 


(3.53) 

(3.54) 


and  ensuing  inequality 


:L  -  o||  :  E  ♦  A  :  ;  -  I  q  •  70  >0 


(3.55) 


where  A  « 

'  3a 

Equation  (3.54)  implies  that  $  is  independent  of  the  temperature 
gradient  76  and  so  is  the  specific  entropy  n. 

We  notice  here  that  we  need  to  know  the  relations  among  L, 

E  and  P  to  obtain  more  restrictions  from  the  inequality  (3.55). 


3.2.1  Caae  I 

We  first  apply  the  kinematical  definition  of  NEMAT-NASSER, 
i.e. ,  from  (3.13) 

t  . 

E-f  FeF  I  dt  (3.56) 

-1  J  ~  ~  1 sym 

PT  *  f1  FPF_1 I  dt  (3.57) 


.  „  ,N  .% 


and  the  inequality 
.pii 


(o  -  p  :Ej+a:Pj+^  :a~—  q’70>_  0 


gives  one  more  restriction 
3  <f> 


a  =  P 


3E 


I 


and  inequality 


!  \PI  +  * :  ““  (T  5’!9  -  0 

where  the  thermodynamic  force  A,  which  is  designed  to  be 
variable  to  the  internal  state  variables,  is  defined  by 
A  = 


3<t> 

3a 


3.2.2  Case.  II 

The  result  (3.15)  of  LEE's  decomposition  gives 

dt 


_  ,t  £e  -e-1 

5n  ■  /  F  f 

o 


|  sym 


■t  ’rP  pP"1 


’  =  F*'  FK  *  -  |  dt 

II  J  ,  'sym 

o 


In  a  similar  way,  we  get 
3  $ 


o  =  P 


3  ET  T 


(3.58) 

(3.59) 

(3.60) 
conjugate 

(3.61) 

(3.62) 

(3.63) 


(3.64) 


(3.65) 


and 


II 


+  A:“'  0 


V8  >  0 


REMARK ;  The  definition  of  plastic  strain  (3.63)  may  lead  to  some 

difficulties, as  noted  by  NEMAT-NASSER  [1979,  1981].  Since  (3.63) 

contains  tensor  Fe  which  represents  the  elastic  deformation,  there 

arises  the  question  as  to  whether  or  not  really  represents 

•  0 

the  plastic  deformation.  Furthermore,  the  definition  F  which  is 

used  to  define  the  elastic  strain  (3.62)  contains  effects  from  the 

intermediate  plastic  configuration  Cp  as  in  Fig.  3.1.  Here  we  assume 

•  • 

that  we  can  vary  the  variables  E  and  P  independently.  When  we 

have  a  deformation  where  significant  couping  between  E  and  P  exists, 
(3.64)  and  (3.65)  may  be  invalid. 3 


3.2.3  CaA<t  III 

From  (3.25) and  (3.26)  we  define  strains  by 
t  •  e  -1 

E-  r  u  u  i  dt 
..III  J  ~  -  | sym 

o 


pm  =  r  upu_1, _ dt 


f  u?u  | 

J  .  -  | sym 


(3.66) 


(3.67) 


then 


'  R  ?m?T  *  ?  !  in?’ 


(3.68) 


f  i~t. 
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Substituting  (3.68)  into  (3.55)  gives 


(RTo  R  -p||-  )  :  Ettt  ♦  RTo  R:P 


3eiii/  *  Jrm  ;  r  r  un 


+  a  :  «-  e  i  0 


(3.69) 


Here  we  rename  R  oR  =  S  and  we  call  S  Dienes  stress  since  it  has 
appeared  in  the  the  work  by  DIENS  [1979). 


By  the  way,  from  (3.69)  we  get 
S  =p-34- 

-  3!m 


(3.70) 


!:!lll  +A:a  -  -  q  •  V  0  >  0 


(3.71) 


REMARK:  Kinematically,  Dienes  stress  S  has  a  special  meaning.  Since 
S  is  defined  as  a  reverse  rotation  of  Cauchy  stress  a,  it  is  the 
stress  in  the  configuration  C  -  in  Fig.  3.2.  The  configuration  Ct~ 
is  designed  to  be  a  rotation-free  configuration.  Therefore,  variables 
U,  S  and  etc.  in  C  -  are  rotation  free  and  we  are  able  to  obtain  a 
materially  objective  stress  rate.  This  issue  has  prompted  many 
discussions  of  appropriate  measures  of  stress  rate  (see,  e.g.  DIENES 
[1979],  NAGTEGAAL  and  DeJONG  [1982],  ATLURI  [1980,  1983]  and  LEE 
and  BERTHEIMER[1983]  ).□ 

From  the  definition 


S  =  R  o  R 


(3.72) 
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3.3  Evaluations  of  Various  Strain  Measures 


Up  to  now,  four  sets  of  system  of  equations  according  to  four 
kinematical  definitions  have  been  proposed.  In  this  section  we 
evaluate  all  of  these  cases  and  choose  a  strain  measure  which  is 
both  physically  sound  and  convenient. 

3.3.1  VeXeAnLna.c.y  o GoveA.ru.ng  Equation* 

Ignoring  thermal  effects  explicitly,  we  again  list  the  following 
equations: 

i)  Field  equations 

p  det  F  *  Pq  (3.80) 

T 

divo+pb=»px,0=0  (3.81) 


ii)  Constitutive  equations 

From  (3.48),  the  plastic  constitutive  equations 


p  -  24. 
3o 


and  the  elastic  constitative  equation  (one  among 


(3.82) 

(3.83) 

(3. 59), (3. 64), 


(3.70)  and  (3.78)) 

o (or  S)  =  p  || 


and  definition 


A 


3 

Ja  ’ 


(3.84) 


(3.85) 


iii)  Kinematical  relations 


The  total  number  of  unknowns  and  the  total  number  of  equations 

are  listed  in  Tables  3.1  and  3.2  for  all  four  cases. 

The  first  and  the  second  cases  give  3  more  equations  than 
unknowns,  i.e. ,  overdetermined  systems  of  equations  are  obtained. 

Generally,  we  cannot  solve  an  over  determined  system  of  equations 
like  these.  So  we  must  find  the  way  to  eliminate  the  extra  3 
equations.  It  is  noted  that  the  3  extra  equations  result  from 
the  introduction  of  the  variable  P  which  may  not  be  a  single  valued 
function.  But  in  the  first  two  decompositions,  we  have  assumed 
that  the  operation  dp/dx  is  valid,  and  this  implies  a  compatibility 
condition  is  fulfilled. 

The  third  and  the  fourth  decompositions  furnish  well  determined 
systems  of  equations.  But  the  fourth  case  Becomes  very  complicated 

since  the  rotation  effects  are  not  explicitly  removed. 

In  the  following  section,  the  actual  differences  of  the 
previously  mentioned  strain  measures  in  the  one  dimensional 
homogeneous  deformation  are  analyzed. 

3.3.2  Calc.alaXA.oM  Elastic.  Stn.aA.ru> 

We  now  consider  a  one-dimensional  case  of  homogeneous 
deformation.  The  deformation  gradients  then  reduce  to  stretch 
tensors  (there  are  no  rotations).  In  this  case,  the  Nemat-Nasser 
decomposition  reduces  to  the  same  representation  as  our  right  and 


SlM 


Table  3.1  Number  of  Unknowns 


Variables  Case  1 


Case  2 

Case  3 

1 

1 

9 

9 

6 

6 

6 

6 

6 

6 

9 

9 

9 

9 

9 

9 

6 

6 

3 

3 

3 

3 

6 

Case  4 


Total 


70 


70 


70 


70 


Table  3.2  Number  of  Equations 


Equations 

Case  1 

Case  2 

Case  3 

Case  4 

3.80 

1 

1 

1 

1 

3.81 

3 

3 

3 

3 

3.82 

6 

6 

6 

6 

3.83 

9 

9 

9 

9 

3.84 

6 

6 

6 

6 

3.85 

9 

9 

9 

9 

3.86 

9 

9 

9 

9 

3.87 

9 

9 

3.88 

9 

3.56 

6 

3.57 

6 

3.89 

9 

3.62 

6 

3.63 

6 

3.90 

9 

3.91 

6 

3.66 

6 

3.67 

6 

3.92 

9 

3.93 

6 

3.76 

6 

3.77 

l - 

6 

Total 


73 


73 


70 


70 
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left  polar  decompositions.  For  simplicity,  we  shall  assume  no 
change  in  elastic  properties  during  plastic  deformation  and  that  the 
material  is  characterized  by  a  bilinear  Cauchy  stress-engineering 
strain  relation  of  the  type  shown  in  Fig.  3.5a. 

Under  these  assumptions,  we  establish  as  possible  configurations 
those  shown  in  Fig.  3.5a.  The  reference,  the  current  and  the  stress 
free  configurations  are  shown  and  were  computed  using  the  definitions 
given  earlier.  An  "elastic"  configuration,  which  is  merely  another 
intermediate  configuration,  is  introduced  for  convenience  for 
calculating  the  Nemat-Nasser  strain  measures. 

Here,  we  consider  a  specimen  whose  initial  length  is  l  and, 
after  a  uniaxial  homogeneous  elastoplastic  deformation  due  to  a  stress 
loading  of  a,  the  specimen  has  elongated  a  length  aXo  as  indicated 
in  Fig.  3.5b.  Then  an  unloading  causes  the  removal  of  the  elastic 
portion  of  the  resulting  total  deformation,  i.e., 


Engineering  strain  £.=(x-X)/X 


(a) 


1)  Reference  configuration  (Cq) 


► 


0  i 


2)Current  configuration  (Ct> 


x=  aX 


3)  Stress  free  configuration  (C  ) 


-  H 

p=  aX-  aoX/(E+c) 


4)  Elastic  configuration  (C  ) 

e 

I  - I 

e=  aoX/ (E+o) 


(b) 

Figure  3.5  (a)  Bilinear  stress-strain  relation. 

(b)  Configurations  in  one  dimensional  deformation. 


Next,  let  us  calculate  Ax,Ap  and  Ae  due  to  a  small  change  in 
o,  so  that  the  increments  (or  rates)  of  various  deformation  gradients 


can  be  obtained: 

x(o  +  Act) 

P(o  +  Ao ) 


„  .  IOOAoX 
aX  +  - E 

„  ,  IOOAoX 

aX  +  - — — 


(a-l-Ao)x 

E+o+Ao 


e  C a  +  Aa  )  = 


(.o+Ao)x 

E+a+Ao 


aE+lQQAo 

E+o-fAa 


Therefore 


Ax  ■  xCa+Aa) 


x(o)_ 


IOOAoX 

E 


(3.98) 


A  p  =  p  (a  +Aa )  -  p  (o  ) 

ao  y  lOOAo  (a -Mo  )x 
E+a  E  E+o+Ao 


(3.99) 


Ae  =»  e(o+Ao)  -  e(o) 
(a+Ao)x  _  ao 

E+a+Aa  E4c 


C3.100) 


We  can  now  compare  the  measures  of  the  two  decomposition 
schemes,  which  are  different  in  this  one  dimensional  deformation. 


i)  Nemat  -  Nasser's  Measures 


3  x 
3X 

D  3  P  010 

s  _  s  CX  —  — 

3  X  E-W 


a 


aE 

E40 


F 


e 


3  e 
3  X 


1+ 


a o 
E+o 


(3.101) 

(3.102) 

(3.103) 


FeAt  -  Fe(o+Ao)-Fe(o) 


3Ae  m  (g+A a)  3x 
3X  °  E+o+Ag  3X 


e  *  p  —  l 

D  At  ■  F  At  F 


(3.104) 


=  AFe  F-1 


(o+Ao)  ..  +  lOOAg  ,  _  g 
E+o+Ag  aE  E+g 

1  ,  100 g  EAg  ,  100g2Ag 

(E+g) (E+o+Ag)  ^  0  aE  aE 


.  100E(Ao)2  .  10Cg(Ag)2^ 

aE  a£ 


As  Ao  *  0  and  A  t  +  0 , 


„e  g  ,  lOOgE  .  100o2^ 

D  •  75?F  (E  +  ~~sr  *  Te  ' 


03. 105) 


oE(E+J )2  CaE2+  100  E°  +  100°2) 


(3.106) 


Similarly, 


^  p  _  3AP  _  gp  lOOAg  (g+Ag )  3  x 

3  x  E-KJ  E  E-K 3  +Ag  3  X 


(3.107) 


DP  -  aE(^a)^  ( 100E2  +  lOOEo  -  aE2) 


(3.108) 


4 


...  ,  ,  (a+Ao)x  ota  „ 

Ae  *  e(a+Ao)  -  e(a)  «  -j- — -4 - -=—  X 

E+a+Ao  E+a 


r  “  3X 


-e  E  +c 
F  =  — 


(3.109) 


AFe  *  FCA  t  =  Fe(o+Aal  -  Fe(o) 

3x(o+Aa)  -e  /  , 

-  !,(.«,)  -  F  <01 


E+a+Aa 

E 


E4a  Aa 
E  *  E 


(3.110) 


DeAt  -  FeAtFe_1 


As  Aa  ♦  0  and  At  ♦  0,  then  we  have 


(3.111) 


.  *  »  -  ^  -  -k  **  .  •  .  “ 


-  **  -  S- 

L«*.JL»»..4L  •*  ?.  J*  — 


ii)  Lee's  Measures 


And  for  D^(or  AF^)  we  have  same  result  (3.107)  [or  (3.108)]. 

Since  the  plastic  strain  measures  in  the  previous  two  cases 
are  of  same  form,  we  now  compare  the  elastic  strain  rates  for  various 
situations.  Table  3.3  shows  the  comparisons  for  various  values  of  a 
with  fixed  ratio  of  o/E  =  .001  and  Table  3.4  shows  the  comparisons 


vv 


Table  3.3  Tabulation  of  proportional  rates  for  o /E 


o 

Proportional 

V  rate 

DeE/a 

.2 

1.4975 

.999 

.4 

1.2477 

.999 

.8 

1.1229 

.999 

1.0 

1.0979 

.999 

1.5 

1.0646 

.999 

2 

1.0480 

.999 

3 

1.0313 

.999 

4 

1.0180 

.999 

5 

1.0180 

.999 

Table  3.4  Tabulation  of  proportional  rates  for  x 


o/E 

Proportional 

rate 

D*E/. 

0 

DSE/  • 

0 

IQ'6 

1.000048 

.999999 

10~5 

1.00048 

.99999 

10'4 

1.0048 

.9999 

io“3 

1.0479 

.999 

io"2 

1.475 

.99 

io'1 

5.37 

.91 

1 

25.25 

.5 

for  various  values  of  cr/E  with  a  fixed  value  of  ct=2.  In  each  Table, 

0  g 

the  values  of  D  E/o  are  tabulated.  In  other  words,  D  E/o=  1  means 
that 

o  =  ED®  (or  ED6)  (3.112) 

3.3.3  Evaluation 

It  appears  to  be  difficult  to  define  deformation,  as  shown  in 
a  convenient  way  that  decompose  into  purely  elastic  and  purely 
plastic  parts  even  in  one  dimensional  deformation  as  shown  in 
Section  3.3.2.  However,  this  is  not  actually  a  necessity  in  a 
meaningful  theory,  provided  the  constitutive  equations  for  each 
are  formulated  properly  (NEMAT-NASSER  [1982],  SIMO  and  ORTIZ 
[1985]).  A  consideration  that  may  be  of  some  importance  is 
whether  or  not  the  particular  decomposition  lends  itself  to  a 

well-determined  system  of  equations  and  to  rotation  effect-free 

strain  measures.  This  proves  to  be  the  case  in  the  new  formulations 

given  here.  A  positive  feature  of  Nemat-Nasser' s  decomposition, 

shared  by  our  own  right-polar  decomposition,  is  that  it  does  lend 

itself  to  experimental  verification  because  of  the  normalization 

by  total  deformation  gradient  Cor  total  stretch  )  in  deformations 
e  •  e  -1  e  #e—l 

(e.g.,  D  *  F  F  i  or  D  =  D  D  i  )  which  may  be  easily  measured 
6  *  -  -  |  sym  1  c'rm 
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Irt 


in  controlled  laboratory  experiments. 

From  now  on,  we  employ  the  following  strain  measures: 


?  •  /  ”Yl|,y»dt 


(3.66) 


!  '  /  PY'lsy  dt 


(3.67). 


CHAPTER  IV 

REPRESENTATIONS  OF  FLOW  POTENTIALS  AND  APPLICATIONS 
4.1  Introduction 

In  Chapter  3,  a  theory  of  elastoplasticity  was  presented 
that  made  use  of  the  notion  of  generalized  stress  potentials  for 
deriving  constitutive  equations  for  the  rate-of-plastic  deformation 
and  the  rate-of -change  of  an  internal  state  variable.  This  theory 
embodies  many  types  of  classical  and  non-classical  plasticity 
theories  and  does  not  require  the  existence  of  yield  functions.  Yield 
functions  can  enter  the  theory,  however,  through  the  definition  of  flow 
potentials,  but  these  potentials  need  not  be  convex  nor  need  they  be 
differentiable  in  the  usual  sense. 

Specific  forms  of  the  free  energy  and  flow  potentials  can  be 
easily  derived  for  cases  in  which  these  are  assumed  to  be  isotropic 
functionals  of  their  arguments.  Such  isotropic  function  representa¬ 
tions  are  by  no  means  unduly  restrictive,  for  we  demonstrate  that 
they  can  be  used  to  represent  some  anisotropic  hardening  rules. 

We  show  that  numerous  special  theories  of  finite  elastoplasticity 
and  viscop lasticity  that  have  been  proposed  in  the  literature  can  be 
deduced  from  such  representations  as  special  cases. 

An  interesting  by-product  of  some  of  these  flow  potential 
representations  arises  because  they  do  not  specifically  involve  a 
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changes  in  elastoplastic  and  viscoplastic  materials  during  various 
programs  of  loading.  As  a  result,  by  defining  a  yield  stress  as  that 
corresponding  to  a  0.2%  residual  strain  (in  a  uniaxial  test)  in  much 
the  same  spirit  as  one  defines  it  in  an  experimental  program,  we 
can  actually  calculate  yield  surfaces  for  such  materials.  Some  of 
these  mathematically  produced  yield  surfaces  are  physically  reasonable, 
others  are  not. 

When  the  representation  of  flow  potential  happens  to  be 
determined  by  a  given  yield  function,  it  is  called  an  associative 
flow  rule;  otherwise,  it  is  called  a  non-associative  flow  rule. 
Classical  theories  of  plasticity  are  generally  characterized  by 
associative  flow  rule  (DRUCKER  1967]  and  PRAGER  [1959],  etc.); 
while  some  authors  (eg.,  ILYUSHIN  [1960]  and  NEMAT-NASSER  [1983]) 
have  shown  cases  in  which  associative  flow  rules  are  inappropriate. 
ILYUSHIN  also  put  forth  arguments  against  the  normality  rule 
which  arose  from  dropping  the  assumption  that  no  change  in  elastic 
properties  takes  place  during  plastic  deformation.  NEMAT-NASSER 
introduced  a  notion  of  a  "workless  strain  rate",  which  could 
invalidate  the  usual  normality  rule  in  calculating  the  rate  of 
stress-work  in  the  isothermal  processes. 

Another  approach  was  advocated  by  HART  [1970,  1976,  1982] 
and  BODNER  and  coworkers  [1972,  1975,  1979,  1983),  who  proposed 
plastic  constitutive  relations  which  involve  a  set  of  internal 
state  variables,  but  not  necessarily  an  explicit  yield  function. 
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Strain  rate  and  rate-of-internal-state  variables  are  independent 
constitutive  variables  in  such  theories. 

A  general  representation  of  flow  potentials  for  initially 
isotropies  materials  with  internal  state  variables  was  proposed 
by  KRATOCHVIL  ]  1978  ]  and  KIM  and  ODEN  [1985a]. 

4.2  Isotropic  Functions 

We  begin  by  reviewing  some  notions  of  isotropic  functions. 

3  3 

2.1:  Let  Q  f  IR  x  ffi  be  an  orthogonal  trans¬ 
formation.  A  function  T  :  U*  1/ ,  where  U  (or  l/)  is  one  of  the  sets 

3  3  3 

(]R,  IR  ,  IR  x  ffi.  ),  is  said  to  be  isotropic  if 

T*(x)  -  T(x*)  V  x  fa  (4.1) 

where  *  is  defined  by 

I)  B*  ■  B,  V  B  (  IR 

ii)  B*  >=  QB  ,  V  B  f  l3 

iii)  B*  -  QAQT  ,  V  B  6  IR 3  x  IR 3 

Next  we  list  the  representation  theorems  for  isotropic  functions 
without  proof.  (For  proofs,  see  GURTIN  [1982],  LEIGH  [  1968]  or 
TRUESDELL  and  NOLL  [1965]). 

Let  G  be  a  subgroup  of  the  symmetric  second-order  tensor  group 
and  SA  be  the  set  of  invariants  of  A,  i.e., 

^ 4  =  {tr  A,  tr  A3,  tr  A3} 


Let 

1(G) 

be  the 

set 

KG)  = 

{IA*  * 

eG> 

Similarly , 

,  we  denote 

IAB  = 

(tr  A, 

tr  B , 

tr  A2 

2  3  3 

,  tr  B  ,  tr  A  ,  tr  B  tr  AB 

tr 

A  B2, 

tr  A2 

B,  tr  A2B2} 

and 

I (G  x 

G)  =  {I 

AB*  A 

x  Bf 

GxG} 

(1) 

Representation 

Theorem  for 

Isotropic  Scalar  Functions 

a)  A  function  <p  :  G  ->•  IR  is  isotropic  if  there  exists  a 
function  <p  :  1(G)  *  IR  such  that 


*(A)  =  0(IA),  7  A  6  G 


(4.2) 


b)  A  function  $  :  G  x  G  +  ffi  is  isotropic  if  there  exists 
a  function  I(G  x  G)  ♦  IR  such  that 


<t>(  A,  B)  =  <J)(  I  Afi)  ,  Va  x  B  (  GxG 


(4.3) 


(2)  Representation  Theorems  for  Isotropic  Tensor  Functions 


a)  A  tensor-valued  function  G  ;  G  *  S  is  isotropic  if 

there  exist  scalar  functions  $q,  0^  ^  (^  )  "*■  ®  such  that 

G(A)  =  0O(  1  A}  i  +  01(  1  A5-  =  ^2(  1  A^2,  ^  A  f  G  (4.4) 

b)  For  isotropic  linear  tensor  functions,  the  function 
S  *  S  is  isotropic  if  there  exist  scalars  y  and  X  such  that 

G(A)  «  2yA  +  A(tr  A) 1  ,  V  A  f  5 


(4.5) 


i  o 


c)  A  tensor-valued  function  G  :  G  x  G  -*■  S  is  isotropic 
if  there  exist  scalar  functions  $  i=  1,2,. ..,8,  :  I  (G  xG)  •"  1R 
such  that 

G (A,  B)  =  1  +  4>  A  +  B  +  <J>  A2  +  <t>  .  B2  +  <P  ,(AB  +  BA) 

+  4>6(A2B  +  BA2)  +  4>?  (AB2  +  B2A)  +  4>g(A2B2  +  B2A2) 

V  (A,  B)€  G  x  G  (4.6) 

It  is  worthwhile  to  mention  the  definition  of  an  isotropic 
material  as  opposed  to  an  isotropic  function  to  perceive  relevancy. 

An  isotropic  material  particle  is  one  such  that  for  some 
reference  configuration  C^,  all  rotations  by  the  response  functions 
of  the  material  are  indiscernable  (cf.  WANG  and  TRUESDELL  [1973]). 

In  other  words,  the  material  is  isotropic  at  particle  p  if  every 
proper  orthogonal  tensor  is  a  symmetry  transformation  for  the 
response  functions  of  the  material  (GURTIN  [1981b],  or  for  a  slightly 
different  version  see  ERINGEN  [1962,  1967]). 

4.3  Representation  of  Flow  Potential  and  Free  Energy 


From  the  theory  of  materials  of  type  N,  we  obtained  the  plastic 
strain-rate  and  the-rate-of-internal  state  variable,  through  relations 
(3.45)  and  (3.48)  respectively  in  the  differential  case,  from  the  flow 


potential.  Next  questions  that  arise  are  how  one  might  obtain  a  spe¬ 
cific  flow  potential  for  a  given  material  and  what  might  be  considered 


as  general  forms  of  flow  potentials  for  certain  classes  of  materials. 

Towards  this  end,  we  shall  first  confine  our  attention  to 
materials  characterized  by  isotropic  functionals.  The  following 
approach  is  similar  to  that  of  KRATOCHVIL  [1978],  although  he  did 
not  make  use  of  the  concept  of  flow  potentials.  We  also  mention 
that  DESAI  et  al. [1983]  suggested  the  use  of  polynominal 
expansions  of  flow  potentials  in  terms  of  invariants  of  its 
arguments. 

From  the  equation  (4.3),  an  isotropic  representation  of  tyfa.A  ) 
is  of  the  form, 

-  2  2  3  3 

v(a,A  )  =  i|>(tr  a,tr  A  ,  tr  a  ,  tr  A  ,  tr  a  ,  tr  A  , 

2  2  2 

tr  crA  ,  tr  aA  ,  tr  a  A  ,  tr  a  A  )  (4.7) 

It  is  occasionally  convenient  in  classical  plasticity  theories 
to  decompose  the  tensors  a,  A  into  spherical  and  deviatoric  parts, 
i .  e .  , 

a  =  si  +  S 
A  =  hi  +  H 

~  ~  (4.8) 

P  =  pl  +  P 

a  -  A 1  +  A 

where  s  =  tr  a/3,  h  =  tr  A  /3,  p  =  tr  P/3  and  A  =  tr  a/3.  Then, 
the  expression  (4.7)  can  be  written  by 


lKo,  A  1  *  ip(.s,  h,  tr  S  ,  tr  H  ,  tr  S  ,  tr  H,  tr  SH, 


-  2  -2  2  -2  2 
tr  S  H  ,  tr  SH  ,  tr  S  IT) 


(4.9) 


The  following  lemma  can  he  easily  proved. 


LEMMA  4.1:  Let  A  and  B  be  second  order  symmetric  tensors  and  be 


decomposed  according  to 


a  *  tr  A/3  ,  A  *  A  -  al 
b  *  tr  B/3  ,  S  =  B  -  bl 


-£«  —  ! 

3A  3  t  • 


3Ctr  I2) 


3(tr  AB)  =  g  ^  3(tr  B4) 


3  (tr  A2B) 


(AB)  +  (BA)  - 


2  Ctr  AB)  1  -  1H 

3  ~~  -  oB 


Hr-CA-1  =  3A2  -  tr  (A2)  1 


3tr(Ol  =  A2  -  i  tr  (l2)  1 
3B  3  ~  - 


3tr(BA2) 


3tr(A2_|2J)  =  (£  g2}  +  (g2  _  l  tr (A  B2)l  =  3 ZlSJLL  1 

aA  ~  ~  ~  -  3  ~  ~  -  A 


When  \  is  differentiable  with  respect  to  both  arguments 


at  (a,  A  )»  by  using  equation  and  Lemma  4.1,  we  can  establish  the 


following  theorem. 


Vs. 
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$ 

*A> 


K;-  •.. 


THEOREM  4.1;  If  the  representation  of  p(0,^  )  is  of  the  form 
of  equation  (4.9)  and  i>  :  w  "*■  1R  is  differentiable  at  (a, A  ),  then 


P  =»  Cftl  +  C,S  +  C-H  +  C,S2  +  C.H2  +  C_  (HS  +  SH) 

U~  1-  2-  3~  4~  5  --  — . 


+C.(§  H2  +  H2^) 

D  «.  „  ^  „ 


(A. 10) 


and 


a  ■  d.l  +  d,H  +  d«S  +  d,H2  +  d.S2  +  dc(HS  +  SH) 
-  1^  2~  3~  4~  5  —  - — 


+d6(HS2  +  S2H) 


(A. 11) 


where  C^’s  and  d^'s  are  scalar  valued  functions  of  the  invariants, 

i-e“  XoA* 

Proof :  The  chain  rule.  Lemma  3.1  and  symmetries  of  S  and  H 
suffice  to  prove  the  theorem.  3 

A  similar  representation  of  the  free  energy  p  of  the  material 
can  also  be  introduced;  e.g. , 


4> CE ,  a  )  »  g>(e,  A  tr  E2,  tr  A2,  tr  e\  tr  t? ,  tr  E  A, 


tr  EA2,  tr  E2A,  tr  E2A2) 


(A. 12) 


where  E  is  given  as  a  sum 


E  =  el  +  E  with  e  =  -j  tr  E 


From  the  relations  (3.61)  and  ( 3. 59) [or  (3. 6A)  or  (3.70)  or  (3.78)],  we 
have 
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A  -  kQl  +  kxA  -  k2E  +  k3A2  +  k4E2  +  k$(EA  +  AE) 

+  kg(E2A  +  AE2)  (4.13) 

0  -  lQl  +  Ijl  +  ^A  +  +  £4 A2  +  £5(AE  +  EA) 

+  £6CA2E  +  EA2)  (4.14) 

As  before,  the  k^'s  and  £^rs  are  scalar  functions  of  the  invariants 
and  density  P.  We  also  note  that,  in  many  cases,  it  is  convenient 
to  introduce  a  decomposition  of  the  free  energy  into  purely  elastic 
and  plastic  parts,  i.e., 

4>CE,  a)  -  <j>e (E)  +  <{>p(a)  (4.15) 

If  (4.15)  holds,  equations  (4.13)  and  (4.14)  reduce  to 

A  -  kQl  +  kx A  +  k3 A2  (4.16) 

;  -  *oi  ♦  *1!  ♦  <4-17> 

REMARK  4,1;  The  inclusion  of  temperature  does  not  result  in  a 
significant  alteration  in  our  procedure.  For  example,  if  we  rewrite 
the  equation  (4.15)  (as  in  KRONER  and  TEODOSIO  .11974])  in  the  form, 
4>(E,  o,  T)  -  <|)e(E,  T)  +  <J>P(a,  T)  (Rl) 

we  would  obtain  same  form  as  (4.13)  and  (4.14).  But  the  coefficients 
become  functions  of  temperature.  ED 

REMARK  4,2;  One  of  the  representations  of  <P  may  be  drawn  from 
linear  elasticity,  i.e.. 
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4>e  =  \  A  (tr  E}2  +  u  tr  (E)2 


(R2) 


where  X  and  y  are  Lame's  constants.  The  elastic  free  energy  (R2) 
results  in  the  relation 


X  (tr  E)  1  +  2y  E 


Since  (S3)  is  invertible,  we  obtain 
E  •»  o  •  r  tr  o  1 

~  fc*  —  £.  — 


which  results  in 


•  1+v  V  v  .  V  , 


(R3) 


(R4) 


(R5) 


where  E  and  v  are  Young's  modulus  and  Poisson  ratio,  respectively, 
and  a  is  an  appropriate  objective  stress  rate  of  Cauchy  stress  o , 
and  E  is  the  elastic  strain  measure  of  (3.66). ZJ 


REMARK  4  .3 ;  It  is  noted  that  in  the  present  theoretical  framework, 
two  representations  of  scalar  functions,  the  flow  potential  and  free 
energy,  are  used  as  compared  to  only  one,  free  energy,  in  some 
thermodynamic  approaches  (e.g.,  KRATOCHVIL  and  DILLON  11969]).  But 
these  approaches  are  equivalent  to  ours  because  one  must  ultimately 
provide  the  evolution  equations  for  the  Internal  state  variables 
and  the  plastic  strain  rate.  These  are  to  be  determined  by  the  flow 
potential  <p  in  our  theory.  □ 


4.4  Examples 


We  now  investigate  how  various  theories  of  plasticity  fit  into 
the  ahove  theory. 

4.4.1  Ca6£  I:  PAa.ndtt-Re.u64  F tout  Ru£e 

The  following  representations  recover  the  Prandtl-Reuss  flow  rule 
which  is  also  obtained  by  defining  the  Levy-Mises  yield  function  as 
a  flow  potential.  Assuming  (4.15)  holds,  set 


-  Kh,  tr  S2) 

(4. 

<t>P  =■  <J>P(A) 

When  the  functionals  in  (4,18)  are 

differentiable  with  their  argu- 

ments,  we  have 

P  -  C.  § 

**  X  — 

(4. 

§  ■  Coi  *  m  co 

(4. 

and 

6  *  koi  ->  h  "  ko 

(4. 

Set  *  A  *  CQ.  Then  (4.19)  reduces  to 

P  =  A  S  (4.: 


which  is  precisely  the  flow  rule  for  the  Levy-Mises  material.  In 
this  case,  h  can  be  interpreted  as  the  bound  in  the  Levy-Mises 
criteria,  i.e., 

f (a)  -  tr  S2  <  a  s  h  »  kn 
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When  h  =  k^(a),  the  representations  (4.18)  describes  a  rate- 
dependent  plastic  material  with  isotropic  hardening. 

4.4.2  Ccu>z  2:  B&ZZ'a  FZouj  Ruiz 
*  -2  2 

If  we  set  A  ■  2  tr  S  /g  ,  where  S  is  constant,  then  we  recover 
Bell's  flow  rule  for  finite  strain  plasticity  (see  BELL  [1983, 
1984]). 


4.4.3  Ca4e  3:  The  PKa.gzn-1-izglzn.  Type.  K-inejncutcc  Ha/idzning 


The  representation 


ip  =  i^(h,  tr  S  ,  tr  H,  tr  S  H) 
<J>P  *  4>PCA,  tr  A2) 


(4.23) 


when  [p  and  are  differentiable,  leads  to  the  following: 


P  =  c:  S  +  c2h 


(4.24) 


a  =  dQl  +  d^H  +  d^S 


(4.25) 


A  =  kQl  +  k  A 


(4.26) 


Then  equation  (4.26)  implies  that 


h  =  k^,  H  =  k^A 


(4.27) 


Introducing  (4.27)  into  (4.24)  and  (4.25)  gives 


P  =  C.S  +  f  A 

M  1  *»* 


(4.28) 


.VN-fo’-  * V •  -  -  *  -  *  •  « 


and 


A  =  f  A  +  d2S 


(4.3 


where  f ^  *  C^k^  and  F2  =  ^1^1  * 

We  now  demonstrate  how  Prager-Ziegler-type  kinematic  hardening 
is  characterized  by  the  representation  (4.23).  The  Prager-Ziegler 
yield  function  may  be  written  in  the  form 

F(!-  sb*  °Y)  =  tr(!  -  ?b)2  -  °Y2  (4-3 

where  is  the  deviatoric  back  stress  which  represents  the  current 
center  of  the  yield  surface.  The  yield  function  (4.31)  results  in 
the  flow  rule , 

P  *  n (S  -  S.)  (4.3: 

«w  -x  D 

where  n  is  a  hardening  parameter.  Next  we  set  ■  -f^  =  d^  *  =  n 

and  =  A.  Then,  equations  (4.28),  (4.29)  and  (4.30)  give  a  set 
of  equations  which  define  a  kinematic  flow  rule  and  an  evolution 
equation  for  back  stress  or  an  internal  state  variable,  respectively 
P  -  ACS  -  A)  (4.3 

and 

A  =  f2  A  +  d2§  (4.3. 

REMARK  4.4:  As  noted  in  equations  (4.24)  -  (4.26)  or  (4.28)  -  (4.30 

p 

the  tensor  valued  functions  D  ,  a  and  A  are  isotropic  even  though 
they  represent  an  anisotropic  hardening  material  (see  equation  (4.6) 
with  appropriate  zero  coefficients).  This  interesting  fact,  that 
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a  non-isotropic  hardening  rule  results  from  an  isotropic  represen¬ 
tation  of  a  flow  potential  functional  or  an  isotropic  form  of  the 

P 

constitutive  equations  for  D  and  A,  suggests  that  references  to 
anisotropic  hardening  may  not  be  proper  terminology.  D 


4.4.4  Caae  4:  Fiouj  Ru£e.  ^oa.  Tjizaca.  McUzsuaZa 


The  Tresca  yield  criteria, 

max  |  a.  -  a.  1-0  <0 

i,  j  1  J  Y- 


(4.35) 


may  be  represented  by  the  non-diff erentiable  functionals 

-2  -2 

i p  =  tjj(h,  tr  S  ,  tr  S  ) 


(4.36) 


4>p  =  4>P(A) 

-2  -3 

The  principal  stresses  are  functions  of  tr  S  and  tr  S  which 
satisfy  the  characteristic  equation 


Oj3  -  Y  tr(S2)ai  -  Y  tr(S3)  *  Q 


(4.37) 


Since  the  Tresca  yield  function  is  convex  but  non-diff erentiable  at 
some  stress  points  (the  comers  of  the  yield  surfaces),  the  classical 
flow  rule  is  characterized  by 


(P,  -a)  £  dip  Co,  A  ) 


(4.38) 


A  "  k0  i 


(4.39) 


where  9ip  ( • ,  • )  is  the  subdifferential. 

9<K<2,A)  =  {(P,  -a)  |  satisfies  the  inequality  (3.46)}. 
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One  can  imagine  generalizations  of  the  Tresca  flow  potential 
to  cases  in  which  is  neither  convex  nor  differentiable.  Consider 
a  hypothetical  material,  the  flow  rule  for  which  is  characterized  by 
a  convex  Tresca  yield  surface  at  a  given  point  in  its  stress  history. 
Suppose  that  low  latent  hardening  prevails.  Then  subsequent  aniso¬ 
tropic  hardening,  after  loadings  in  two  distinct  direction  as  in 
Fig.  4.1,  could  conceivably  produce  a  non-convex  yield  surface  of 
the  type  shown  in  Fig.  4.1.  While  some  may  argue  that  such  hardening 
patterns  are  rare  or  non-existent  in  real  materials,  they  are  never¬ 
theless  acceptable  in  the  general  theory  outlined  here.  As  an 
example  of  non-convex  yield  surface,  WILKINS  et  al  11980]  derived  a 
star  shaped  yield  surface  such  as  that  in  Fig.  4.2  from  experiments 
on  metal  specimens.  D 

REMARK  4.5 :  We  also  remark  that  there  is  frequently  some  ambiguity 
in  the  literature  as  to  exactly  what  is  meant  by  the  yield  surface  of 
a  given  material.  BELL  11983,  1984]  uses  two  yield  surfaces,  the 
inner  surface  and  outer  surface,  and  HILL  J1979J  also  has  distin¬ 
guished  between  surfaces  corresponding  to  an  elastic  limit  and  a 
plastic  limit.  In  subsequent  examples, we  use  a  single  yield  limit 
(or  surface),  as  the  elastic  limit,  for  simplicity . G 

Now  let  us  investigate  a  proposed  new  model  of  plasticity  (or 
viscoplasticity)  which  need  not  employ  the  concept  of  a  yield 


surface . 


4.4.5  CoAe.  5:  BocLma.  eX  at.  ftou)  RuZe. 


4.4. 5.1  JiotAopZc  HctA.de.nLng 


BQDNER  and  PARTOM  (1972,  1975]  proposed  a  flow  rule  which  is 
drawn  from  the  classical  Prandtl-Reuss  flow  rule  and  improved  by 
phenomenological  observations,  but  which  does  not  involve  specifi¬ 
cation  of  a  yield  surface.  This  theory  is  offered  as  a  way  for 
characterizing  isotropic  hardening  of  certain  materials. 

For  the  Bodner/Partom  theory,  the  following  form  of  potentials 
can  be  introduced: 


ip  = 


i=0  (2ni-l)i!  J2(ni-1,) 


(4. 


<J>P  =-11^-  hQ(h1-h0)  exp  (-mX/h0)/m 


(4. 


Here 


J2  -  j tr  !2 


B2  =  4  h2  C— ) 
J  n 


l/n 


and  Dq,  h^,  h^  and  n  are  material  constants. 

From  the  expressions  (4.40),  (4.41)  and  the  generalized 
normality  rule  (3.45)  (in  this  case,  equation  (3.48)),  we  obtain 
the  following  Bodner/Partom  flow  rule. 
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* 

1 

I 

I 

■ 


with 


P  =  (Dq  exp[-B2/J2)n^)S/J2l5 

* 


(4.42) 


fa  =  +  (faQ  -  fa1)  exp  ( — mA/ VIq) J 


It  is  noted  that  as  n  +  «,  equation  (4.421  models  non-work 
hardening  materials  (i.e.,  tfae  elastic  perfect  plastic  material). 


4.4. 5.2  Non-lAo&iopi.c  HaA.de.ning 

STOUFFER  and  BODNER  X 197 9]  extended  the  previous  hardening  model 
to  the  case  of  non-isotropic  materials.  A  plasticity  theory  for  non¬ 
isotropic  hardening  can  he  constructed  by  introducing  implicit 
forms  of  potentials  which  are  different  from  the  previous  isotropic 
representations,  I.e., 

^  -  ip(o.  A,  ®21  (4.43) 

aP  =  a(a^, 

where  £  C^(1R^  x  IR^  x  IR^l  and  o  £  C^(IR^  x  IR  ^),  and  we  employ 
two  internal  state  variables,  and  a  .  As  one  example,  suppose 
p 

that  D  is  given  by  a  generalized  Prandtl-Reuss  flow  rule 


or 


!  •  f 


P  ij  =  KijU  J 


-1  =  '3o7  *  po  i  ?  h  =  W 

*2  ‘It’  pi(51i  *  V!!1?: 


( 4 . 44  ) 


(4.45) 


« 

i 


The  symmetry  of  D  and  O  and  the  reciprocity  condition  of  HILL  1 1979] 
result  in  the  conditions 


^ijk£  ik£  ^ij£k  ^k-tijj 


(4.46) 


If,  in  addition,  we  assume  that  the  hydrostatic  component  of  the  stress 
does  no  work,  then  equation  (4.44)  becomes 


P-KCAilSorDij  -KijR  Sk 


(4.47) 


(4.47) 


Due  to  conditions  (61),  equation  (62)  can  be  rewritten  in  the 


form 


d  =  6  T 


(4.48) 


where 


Next,  a  principal  hardness  direction  can  be  defined  such  that 

d^  =  g  T  (a  is  not  summed)  (4.49) 

a  a  a 

so  that  we  can  rewrite  equation  (4.45)^  in  a  form  similar  to  (4.49): 

A  =  P.  (h)  +  P.(h)ct 
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Based  on  experimental  observations,  STOUFFER  and  BODNER  [1979] 
provide  the  following  empirical  relations  for  each  term  in  (4.49) 
and  (4.50); 


h  '  °o  “p  =£  <A„2  h” 

A1  =  hn  +  9  f  *»(T)  dl  +  (l-q)a  f  Yni(T)h(T)dx 
a  0  JQ  %  “ 

h  «=  mChj-hg)  expl-m  X/hQ]X/h0 


(4.51) 

(4.52) 

(4.53) 


with 


(4.54) 


*  ■  ?  !  5P '  E  \  dl  <4-55> 

a 

Here  the  coefficients  DQ,  n,  h^,  h^,  m  and  q  (the  latter  specifying 
the  degree  of  isotropic  hardening}  are  material  constants,  Ya 
represents  the  cosine  of  the  angle  between  the  loading  direction  and 
the  principal  hardness  direction,  and  ||  •  |[  denotes  the  Euclidean  norn 


4. 4. 5. 3  CaZc.uZaXA.on.  orf  VieXd.  Sun-^acex 

An  interesting  by-product  of  the  plasticity  theories  discussed 
here  is  that  one  can  calculate  an  equivalent  yield  surface  by  setting 
an  arbitrary  permanent  strain  level  as  an  indication  of  yielding.  In 
fact,  this  rather  arbitrary  selection  of  a  measure  of  yielding  is  no 
different  than  what  is  done  experimentally  to  determine  yield  stresses 


for  most  metals. 
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In  the  present  cases  we  suppose  that  we  are  supplied  with  a 
specimen  of  material  characterized  by  a  collection  of  constitutive 
equations  of  the  Bodner-type  Ie.g.  (4.51)  -  (4.55))  and  that  an 
essentially  uniform  state  of  stress  is  developed  during  a  laboratory 
test  (e.g.,  a  uniaxial  tension  test  of  a  cylindrical  bar  of  the 
material) . 

We  shall  assume  that  yielding  of  the  material  occurs  whenever 
a  .2  percent  offset  permanent  strain  is  experienced.  We  shall  thereby 
trace  out  an  effective  yield  surface  for  this  material  by  determining 

the  yield  stress  in  various  directions  and  connecting  these  stress 
points  in  the  ir-plane 

A  numerical  procedure  for  performing  this  construction  is 
outlined  as  follows: 

A)  Numerical  Scheme 

An  explicit  method  for  integration  of  rate  type- variables  is 
used  for  convenience.  The  following  algorithm  is  used: 

(1)  specify  as  an  increment  of  stress, 

Tn  =  n  T  /N  (4.56] 

where  N  is  total  number  of  Increments  and  T  is  the  maximum  applied 

-max  rr 

load. 


(21  Calculate  the  increment  of  the  hardness  variable  h, 
Ahn  =  m(hQ-h1)  exp[-m  An  1 ] AAn  */hQ 


(4.57 


* 


AAn_1  =  I  TI)"1APn-1  and  X11”1  =  Z  AAk 
i  1  1  k=l 


C3)  Calculate  ^ 


a 


A"  =  h  +  q  E  Ahfc  +  CYq)  slgnlAc^]  I  Ahnyn 
u  k-1  1  fc=l  a 


(A. 59) 


(A)  Calculate  the  plastic  'moduli 


(A. 60) 


(.51  Calculate  the  plastic  strain  Increment 


*  _n  _n 

AP .  *  6  T 
i  a  a 


(A. 61) 


Next  we  proceed  to  the  calculation  of  the  yield  stress.  During 
this  calculation,  the  hardness  variables  A^  are  assumed  to  be 
constant  since  up  to  .2  percent  plastic  strain  there  should  be  no 
significant  change  in  A^. 

(6)  Increase  the  stress  in  various  directions  until 
.2  percent  of  an  appropriate  measure  of  total  strain  ||pj|  is  attained 


e^-^(Aa/3j;i]J2 


(A. 62) 


(7)  Project  the  yield  stress  into  the  ^-plane  as  follows: 

Cal  First,  take  stress  vectors  which  lie  in  the  - 
c?2  plane,  since  we  change  the  direction  of  stress  vector  in  the  - 
02  plane  for  convenience  In  step  (6).  Then  calculate  the  magnitude 


of  that  vector  and  the  angle.-i»etween  the  vector  and  the  axis 
(see  Fig.  4.31: 

a  =  |  a  |  =  Vax2  +  a2 

9  *  Arc  tan  (a^/a^l 


Cbl  Project  a  into  the  intersection  line  and  the 
line  perpendicular  to  the  intersection  line  which  lies  “  ”  <?2 

plane, 

os  -  |  a  |  sin  (J  +  01 
a  *  | a |  cos  Cx  +  91 

(cl  Project  0g  into  ir-plane  and  compute  the  magnitude 
of  the  projection  0  of  a  , 

TT  S 


and 


°s  C0S  \ 


a 


c 


where  <p  ■  cos  ^  (1/yTT 

(d)  Compute  the  angle  between  the  line  perpendicular 


to  the  intersection  line  in  the  u-plane, 


Figure  4.3  Projection  of  the  stress  vector  into  the  -n-plane. 
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<f>  f  -  cos  1  (4>_  /O  ) 

I  IT  1  IT 

Throughout  the  steps  (a)  -  (d),  the  yield  stress  can  be  placed  in 
the  it  -plane. 

B)  Numerical  Experiments 

First,  we  consider  a  specimen  of  titanium  RMI-50A.  The 
material  constants  of  this  specimen  listed  in  BODNER  and  PARTOM 
[1975]  are: 

m  =  100 
n  =  1 

Dq  =  104  sec  " 1 

hQ  *  1150  N/mm2  («  167.9  Kpsi) 
h^  =  1450  N/mm2  (*  211.7  Kpsi) 

The  elastic  constants  are  taken  to  be 
K  =  123,000  N/mm2 
G  =■  44,000  N/mm2 

To  check  our  algorithm,  the  stress-strain  curve  produced  in 

2 

a  uniaxial  test  with  the  maximum  loading  T  =  400  N/mm  was 

max 

calculated  and  the  results  are  shown  in  Fig.  4.4. 

Next  a  yield  surface  is  created  after  a  uniaxial  loading 
history  for  the  case  q  ■  1  and  q  *0  as  indicated  in  Fig.  4.5 
and  4.6,  respectively.  In  Fig.  4.5,  we  recognize  that  the  Bodner 
et  al  theory  is  equivalent  to  the  Levy-Mises  yield  criterion  in  the 
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case  of  q  *  1.  Quite  different  results  are  attained  in  the 
hardening  case  illustrated  in  Fig.  A. 6;  not  only  do  corners  appear, 
but  a  marked  lack  of  convexity  is  seen  to  appear  as  the  yield  surface 
assumes  a  "key-hole"  shape. 

4. A. 6  Ccu>&  I/I:  MocU.4-ie.cl  BodneA.  and  Pa/Uom  Theory 

In  the  previous  section,  we  analyzed  the  constitutive  equations 
by  Bodner  and  his  coworkers.  In  their  theory,  the  variable  X 
conjugate  to  the  hardness  variable  h  needed  to  be  defined  but  in 
our  theory  of  materials  of  type  N  this  variable  can  be  derived 
from  the  flow  potential  ip.  Here  we  will  furnish  the  theory  which 
follows  the  theory  of  materials  of  type  N,  that  implies  thermo- 
mechanical  validity,  and  we  will  use  these  constitutive  equations 
for  the  finite  element  computations  in  the  next  chapter. 


A. A. 6.1  Equation* 

We  first  introduce  a  pair  of  potentials  as  in  the  previous 


(A. 63) 

(A.6A) 


1  n+ 1 

Where  B  =  ■—.(—  )  and  X,  p  are  Lame's  constants  and  D  ,  h  ,  h,  and  n 
a  n  o  o  1 

are  material  constants. 

Applying  (3.48)  and  results  (3.61)  and  (3.70)  to  the  potentials, 
we  obtain  the  following  elastoplastic  constitutive  relations. 


9<|>  pA  p 

' 0  a?  *  p-  <trE)1  2uE 

0  0 

(4.63) 

3  4* 

=  “  7x  =  hl  +  (ho‘hl)  exP(-mX) 

(4.66) 

3i|»  „  1  ,  Bh2", 

-  9S  Do  h/J2  exP  <-  ,  n  >  ? 

J2 

(4.67) 

.  _  ii  .  I  s’-  p 

3h  h  .  . 

(4.68) 

where  S'  =  S  -  -j  tr(S) 


From  the  thermodynamic  restriction,  we  obtained  the  inequality 
(3.71).  We  need  to  check  whether  the  given  constitutive  relations 
satisfy  this  inequality.  Otherwise,  we  have  to  impose  this  inequality 
to  get  a  correct  set  of  solutions. 


S  : 


[D°  h/Jl  SXP  (" 


Bh 


2n 


-)]  S' 


=  [positive  value]  (2J2) 


_3$ 

3a 


a 


3<t>  • 

Tx  :  A 


=  hX  =  S'  :  P  >  0 


Thus,  -  (1/9)  q • V9  _>  0,  as  generally  required. 


4. 4. 6. 2  DzXz/un^.natAon  06  Coe.66j.CA.zntA  (Ja om  £ x.peAAmzntal  Rzau ItA 
We  now  have  a  complete  set  of  constitutive  equations  for 
elastoplasticity .  A  natural  question  that  arises  here  is  how  to 
determine  the  material  constants  experimentally. 

To  determine  the  elastic  constants  X  and  u  (or  Young's  modulus 
E  and  shear  modulus  G),  we  need  two  experiments,  i.e.,  tension 
test  and  torsion  test  in  the  infinitesimal  elastic  strain  range. 

For  the  plastic  parts,  there  are  5  constants  to  be  determined, 
i.e.,  Dq,  hQ,  hj,  m  and  n.  To  determine  these  constants  we  need 
at  least  two  uniaxial  tension  (or  compression)  tests  with  different 
strain  rates  since  the  constitutive  equations  have  a  capability 
to  model  rate  sensitivity.  Next,  let  us  examine  the  role  of 
individual  coefficients. 

A)  Dq  acts  as  a  scaling  factor  to  the  plastic  strain  rate. 

It  also  effects  the  yield  stress  in  a  classical  sense. 

B)  hQ  and  h^  are  minimum  and  maximum  values  of  hardness  variable. 

h  will  mainly  determine  the  yield  stress  and  h.  will  limit  the  maximum 
o  i 

achievable  stress  (actually  J2). 

C)  m  adjusts  the  growth  rate  of  the  hardness  variable  h  with 

respec t  to  X  . 

D)  n  determines  the  hardening  tendency  of  the  material.  We  may 
obtain  a  elastic-perfect-plasticity  as  n  -  °°  .  Usually  n  =  1 

is  chosen  for  the  hardening  material  since  the  detailed  approximation 
can  also  be  done  through  the  previous  coefficients. 


Next  we  explain  how  to  determine  those  material  constants 
from  experimental  data.  The  stress  strain  curves  in  Fig.  4.7  were 
taken  from  the  experiments  done  by  Instron  1125  test  machine  with 
pure  lead  (more  than  99%)  specimens  which  have  25  mm  gauge  length 
with  cross-sections  size  3.8  mm  x  6.5  mm. 

An  Instron  extensometer  was  used  for  strain  (in  fact,  the 
total  displacement  of  gauge  length)  measurement  and  the  load  was 
recorded  as  a  function  of  the  strain.  The  crosshead  velocities 
applied  were  0.5,  5  and  200  mm/min.  which,  for  an  effective  over¬ 
all  gauge  length  of  25  mm,  correspond,  respectively,  to  the  constant 
engineering  strain  rate  of  3.3  x  10  4,  3.3  x  10  3,  and  1.33  x  10  l, 
respectively.  The  experimental  curves  at  the  lowest  and  the  fastest 
rates  are  chosen  to  determine  the  material  constants. 

First  we  choose  n  -  1  as  a  common  choice  for  hardening  material. 
Let  a  uniaxial  stress  in  the  experiment  be  o,  i.e. 

°x  -  a  ,°y  ■  0  ,  °z  =  0  (4.69) 


and 


=  20/3 


»  _  5. 


»  a. 


(4.70) 


III  ♦  I  ♦  I)  a2  -  £2 

2l9  9  9J  0  3 


Thus  the  constitutive  equations  (4.65)  -  (4.68)  in  this  state  of 
stress  become  simplified  as  follows: 


1.32  x  10_1 

1.27  x  10_1 

10.58  N/mm2 

3.3  x  10‘4 

3.2  x  10'4 

8.73  N/mm^ 

Then  from  this  table  and  (4.72),  we  have 

Dq  -  1029400  N/mm2  sec 

h  =  26.73  N/mm2 
o 


To  determine  h^  and  m,  we  need  to  calculate  h.  Taking  logarithm 
of  (4.73)  gives 


In  ?1  «  In 


2^3D 

3 


-  £n  h  - 


(4.75) 


By  using  the  Newton  Raphson  scheme,  we  obtain  the  variation  of  h 
versus  strain  as  in  Fig.  4.8. 

With  the  computed  value  h,  we  are  ready  to  determine  h^  and  m 
First,  take  time  derivative  to  (4.72). 

h  =  -m(h  -h^)  exp  (-mX)A  (4.76) 


And  then,  from  (4.72)  and  (4.76)  we  obtain 


ENGINEERING  STRAIN.  % 


Figure  4.8  Computed  variation  of  the  internal  state  variabl 
from  experimental  data 


ENGINEERING  STRAIN.  % 


Figure  4.9  Experimental  and  calculated  stress  and  strain  curve 


CHAPTER  V 


FINITE  ELEMENT  APPROXIMATION 

In  this  chapter,  we  solve  a  boundary  value  problem  with  the 
constitutive  equations  introduced  in  the  Section  4.4.6.  We  limit 
ourselves  to  quasi-static  problems  involving  no  body  forces  and 
no  temperature  dependence.  That  means  we  do  not  solve  the  energy 
equation  directly  and  that  the  equation  of  momentum  balance  in 
the  reference  configuration  will  be  of  the  form 

V  T  -  0  in  C  (5.1) 

x  .  .  o 

where  T  is  the  first  Piola  Kirchhoff  stress  and 

7  M  i- 
x  3x  . 

Early  accounts  of  finite  element  models  of  elastoplasticity 
were  reported  by  ODEN  and  KUBITZA  [1967],  MARCAL  and  KING  [1967], 
ARGYRIS  [1967],  and  ZIENKIEWICZ  et  al  [1969],  and  finite  defor¬ 
mations  were  treated  by  ODEN  [1968,  1970,  1972]  and  others.  A 
variety  of  different  formulations  of  the  large  deformation  problem 
have  been  explored,  and  we  mention  the  incremental  Lagrangian 
formulation  of  HIBBIT,  MARCAL  and  RICE  [1970],  the  updated 
Lagrangian  scheme, of,  e.g.,  MCMEEKING  and  RICE  [1975],  and  various 
related  schemes  proposed  by  NEEDLEMAN  [1972],  LARSEN  and  POPOV 
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[1974],  BATHE,  RAMM  and  WILSON  [1975],  ARGYRIS  et  al.  [1977,  1978], 
CESSCOTTO,  FREY  and  FOUNDER  [1979]  and  ATLURI  [1980].  Numerous 
applications  of  finite  element  methods  to  metal  forming  problems 
can  be  found  in  the  literature  and  we  mention  as  examples  those  in 
the  papers  of  LEE,  MALLET  and  YANG  [1977],  NAGTEGAAL  and  DeJONG 
[1980],  KEY,  KRIEG  and  BATHE  [1979],  ARGYRIS  and  DOLTINIS  [1979, 

1980],  TAYLOR  and  BECKER  [1983],  and  KIKUCHI  and  CHENG  [1983].  A 
good  survey  of  current  theories  and  numerical  methods  for  finite- 
deformation  plasticity  can  be  found  in  the  volume  edited  by 
NEMET-NASSER,  ASARO,  and  HEGEMIER,  [1984]  and  in  the  proceedings 
edited  by  WILLAM  [1984];  see  also  the  recent  work  of  SIMO  and 
MARSDEN  [1984],  SIMO  and  ORTIZ  [1985],  and  KIM  and  ODEN  [1985] 
and  the  references  therein. 

5.1  Formulations 

Here  we  adopt  the  incremental  total  Lagrangian  formulation, 
where  the  reference  configuration  is  always  the  initial  configuration. 
Let  q  denote  a  test  function  which  belongs  to  a  set  of  V  admissible 
displacement  increments.  Then  a  weak  form  of  equation  (5.1) 
can  be  written 

f  T  :  Vq  dx  -  f  TN-q  dS,  f  v  (5-2) 

J  no .  J3$20"  ' 

where  N  is  the  outward  normal  unit  vector  in  the  reference  configu¬ 


ration. 


We  assume  that  the  equilibrium  state  is  achieved  at  n-1  th 


incremental  step,  i.e., 

JV»-i;75  ?  •  ? ds-  V?<fv 


(5.3 


Then  at  the  nth  step. 


n 

r  *  l  AT.  -  T  ,  +  AT 
.n  .  ,  .i  .n-1  .n 
i*l 


t  =  T  N  »  t  .  +  At 
_n  .n.  _n-l  _n 


and 


■  4)  ds-  Vqfv  (5.4 


From  the  definition  of  the  first  Piola-Kirchhof f  stress  and 
equations  (3.68)  and  (3.72), 

AT  »  A(JRSU_1) 

»  A(A(trE)F_T  +  2PREU*1) 

-  A(trAE)F_T  +  2 yRAE  U_1 

+  A(trE)AF~T  +2uAREU-1  +  2uREAU_1  (5.5) 

By  applying  the  decomposition  relation  (3.68)  and  the  form  (5.5), 
equation  (5.4)  becomes 
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of  (5.7)  are  obtained  by  taking  &P,  AR  and  AU  as  functions  of 
U  and  transferring  these  to  the  right-hand  side,  i.e.. 


.  .  o AP  i-1  .  3AR  i-1  . 

r  *  tttt  AU1,  AR  =  — 77  AU1 ,  etc. 

_n  3AU  n  ,n  .n  3AU  _n 

n 


Note  that  AP^  ,  etc.  go  to  zero  when  the  iteration  converges 

since  P  ,  R  ,  etc.  are  defined  in  the  following  fashion: 

_n  _  n 


P  =  P  ,+ZAP  ,  etc ... . 
,n-l  .  _n 
l 


Adding  (5.3)  to  (5.7)  gives 


J\  [ Xtr ( At  D1)  ♦  2 u  At  D1  ]  (  f'V"1  :  yq  dx 
,n  .n  _  n 


/  [X(trAP)F~T  +  2p(R  AP  U_1)]1_1:Vq  dx 

“o  -  '  '  '  '  n 


/  [  X  ( trE)  AF”"T  +  2u(  AR  E  U_1  +  R  E  AU_1)]1_1:Vq  dx 

“o  '  '  -  -  -  '  -  '  n  ”  ' 


V-  q  V 


J' 


FI  F2 


F3  F4 


and  !V  »c 


(5.11) 


F1Vl  +  F2^,2  -  6A 


F3*8,l  +  F4*B,2  “  6B 


AtD1  = 
~n 


Y  { 3 B  Ug  +  BA  U2} 


Y  <SA  u2  +  SB  Ug } 


(5.12) 


To  reduce  bookkeeping,  we  note  that  the  terms  AFd  *  and  (AU  *) 

-n  -  n 

should  tend  to  zero  when  convergence  in  each  iteration  procedure 
is  achieved.  It  thus  would  allow  us  to  delete  the  second  integral 
on  the  right-hand  side  of  (5.8).  Then  the  discretized  equations 
of  (5.8)  at  each  node  in  nth  element  would  assume  the  form. 


,  ,  ,  fcuO  r(Fl*  ,BA-i(F3*  ,BB*F1*  ,BB)}u' 

l,1.,2)  /jWxTRd)  *  2„  0-‘  2  “•1  “-2  f 

0  bBJ  L(F2»0i16A  »2<F4*0i1SB*F2»  28B))us 


.  (i(F3*  .BA  *  FI,  BA)  ♦  F3»  ,BB)uf 
2  <3,1  a, 2  a, 2  3 


♦  { 4(  F4^)  BA  +  F2<p  ,6A)  +  F4$  ,6B}u.2^ 

L  &  *  ■!•  Cl  *  ^  Cl » p 
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{q1 >q2}  \  J  |  x(pi  ♦  P2) 


(tx  +  Atx) 


a  a]  fell*  .  ♦  C21*  ,1 
.. :  “.I  o»2 


(t2  +  At2) 


aB  C12$  .  +  C22 <t>  ,  I 

o,l  o ,  2  |  / 


RS U 1 4>  ,  ♦  RSU3*  , 

a,l  a, 2 

M  RSU2*  ,  +  RSU4$ 

N  a,l  a, 2 


dx  \  , 


(5.13) 


where  ^  is  an  nth  element  and  is  the  side  of  the  element  on 
which  the  traction  is  prescribed.  Summations  on  a  and  B  are 
implied.  Also, 


.  .  f  R1  R3l  .  ,  V?1  P3l  .  .  f U1  U3 

[R]"-1  =  I.UPli"1-  .  tujJ-1-! 

'  |_-R3  R2  I  [_P3  P2j  '  [_U3  U2j 


2u 

Det  U 


RP1 *U2-RP2  »U3-RP1 *U3+RP2  «U1 

REP.U2-RP4.U3-RP3.U3+RP4.U1 


,  [RP]  =  [R]  [AP] 


I1'1  -P 


SI  S3 
S3  S2 


,  ( RSU  ]  * 

RSU1 

RSU  2 

RSU3 

RSU4 

[R]  [S][U_1 ] 
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and 


aA  =  FI#  +  F3 d>_  - 
a,l  a, 2 


aB  =  F2d>  +  F4#  . 

°,1  a*2 


Again,  by  setting 

aC  =  F3#  .  >  FI*  , 

Cl  j  1  Cl  y  2 


a.D  =  F4#  ,  +  F2#  - 

a,l  a, 2 


equation  (5.13)  can  be  reduced  to  the  following  set  of  simultaneous 
equations  for  each  element: 


6r  62  ’  ‘ ' '  6nn 


K12_ 

< 

"  1  1 

u6 

» 

2 

_2 

L*21 

K22 

—i 

ufi 

J 

(5.14) 


for  every  a  -  .  °2 » • • • °^N- 


Here  the  a^'s  and  6^'s  are  global  node  numbers  in  an  element  and 
NN  is  the  number  of  nodes  in  each  element. 

In  global  form, 

N 


I  .  I 


8. ,  6. . . .  6 


1 

K12! 

L21 

K  ?2  1 
J 

C  , 

r  ^ 

1  1 

_  i 

U  8 

4 

>  < 

) 

1  2 

2 

u  a 

'  F 

(5.15) 


for  every  a 

where  N  is  the  number  of  elements,  ELEM  is  an  element  number  and 

a8K,.  *  f  .,{A*aA*8A  +  2p,Fl»<f>  ,BA  +  p*aC*BB}  dx 

11  J _N  a,l 

“o 

aBK.»  *  f  „{A*aA*8B  +  2u * F3* _BB  +  u*ctC*8A}  dx 

12  ■>  0N  a, 2 

0 

aBK.,  ■  f  „{A*aB‘BA  ♦  2p * F2*<t>  ,BA  +  p*ctD*BB}  dx 

21  JaN  ot,l 

0 

a6K._  =  f  „{A*aB*8B  +  2u*F4*$  -BB  +  p«aD*BA}  dx 

22  JjjN  a, 2 

0 

F»  -/0N(<A<P1  *  P2)F1  *  C11  -  RSU1>*a,l  *  (><P1  *  P2)F3 
0 

.  C  -  K8U3*  Hx  *  /  t[  ds 

dQa  n 

Fa  =  /  N  {(X(P!+P2>F2  +  Ci2  "  RSU2Ha,l  +  (X(P1  +P2)FA 

no 

+  C__  -  RSU4)$  }  dx  +  f  tj  <|>  ds 

22  0,2  3fi  2n  a 

a 

REMARK : 

i)  In  the  first  iteration  at  each  incremental  step,  we  impose 
the  incremental  essential  boundary  conditions:  but  after  the  first 
iteration,  we  must  impose  the  zero  values. 

ii)  It  is  interesting  to  note  that  (5.15)  reduces  to  equation 
for  stiffnesses  in  linear  infinitesimal  elasticity  upon  appropriate 
specializaton.  In  the  case  of  infinitesimal  deformation, 
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F  =  1,  R  =  1,  U  =  1 


aA  =  $  ,  ,  aB 

®  *  1 


,2  ’  aC  =  *0,2  ’  aD 


6A=  *8,1  ,  68  =  *8,2 

Therefore,  as  expected, 


K11  ■  /2SI<X*  2ll>*«, 1*8,1  *  "‘a, 2*8, 21  d? 

’*12  ‘  /dIX*a,  1*6,1  *U*a,l*B,2l  dl! 

>K21  -/aN[Ha,2*S,l  '“♦a, 1*8,2'  d? 

'K22  ■  /„N[<>  *  2B)*a,2*8,2  *  U*a,l*B,l'  “? 


Note  that  the  assembled  matrix  then  becomes  symmetric. 

In  each  iteration,  the  following  constitutive  routine  has  to 
be  solved: 

h^  *  h1  ♦  (hQ  -  h^  exp  ( -mfv ^  *) 

i  •  i  i-1 

AP^  =  A  t  P(h~j ,  0  L) 

-  n  n_  n  ,  n 

-i  i 

AE1  =  AU  U"1  |  -  APn 

-  n  ,  .  sym  n  ,  l 


AS1  =  -  [XCtrAE1)!  +  2UAE1] 

- n  jti  -n  -  -n 

AA  »  i-  S  :  AP1 
n  ,1  -n  -n 
h 

n 

In  the  actual  computation,  we  subdivide  the  incremented  solution 
Au  U  ^|Sym  by  a  prescribed  number  and  proceed  by  using  the  previous 
forward  Euler  method. 

Since  the  previous  constitutive  equations  assume  incompres¬ 
sibility  in  the  plastic  deformation,  care  must  be  taken  in  choosing 
a  stable  element  approximation  (see  NAGTEGALL,  PARKS  and  RICE 
[1974]  and  the  detailed  stability  analysis  of  ODEN  et  al  [1982, 
1984]).  Here  rectangular  elements  which  consist  of  four  3-node 
triangles;  so-called,  four  constant  strain  triangles  (4CST  element) 
issued.  A  mathematical  analysis  of  this  4CST  element  can  be  found 
in  KIKUCHI ,  ODEN,  and  SONG  [1982]  and  KIKUCHI  [1983]. 

5.3  Numerical  Examples 

In  this  section  several  example  problems  are  solved  to  verify 
the  algorithm  described  above. 

The  proposed  set  of  constitutive  equations  (equations  (4.65), 


(4.66),  (4.67)  and  (4.68)  contains  7  constants  to  be  determined. 


from  experimental  data.  The  two  elastic  constants  can  be  determined 
by  standard  procedures, but  the  plastic  constants  require  at  least 
two  uniaxial  tests  at  different  strain  rates. 

Following  Bodner  and  Partom,  we  consider  estimated  constants 
on  Titinium  for  our  constitutive  equations,  which  are  similar 
to  theirs.  The  material  parameters  are  as  follows: 

i)  Elastic  Constants 

X  =  93667  N/mm2  (1  N/mm2  =  106  Pa) 

V  *  44000  N/mm2 

ii)  Plastic  Constants 
n  *  1 

m  *  50 

D  -  1.35  x  107  sec'1 
h  =1150  N/mm2 
h  =  1450  N/mm2 

As  a  first  computational  example  the  homogeneous  plane  strain 
elongaton  of  a  test  block  is  computed.  Results  are  shown  in 
Fig.  5.1.  The  computed  results  reflect  up  to  20  percent  engineering 
strain  and  show  strong  sensitivity  to  strain  rates.  Note  that 
slower  loading  results  lower  yield  stress. 

A  result  of  loading,  unloading  and  reloading  is  shown  in  Fig. 
5.2  with  a  strain  rate  of  1.5  x  10  3  per  second.  The  so-called 
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Stress-strain  curve  with  loading-unloading 
at  a  strain  rate  »  1.5  x  10~3. 
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ratchet  effect  at  yield  point  in  the  reloading  process  is  not 
apparent  because  of  a  large  amount  of  deformation. 

In  the  third  example,  the  variation  of  the  internal  state  vari¬ 
able  with  strain  is  shown  in  Fig.  5.3.  In  this  case,  the  internal 
state  variable  represents  hardness  of  the  material.  It  is  seen  to 
vary  with  strain  in  a  way  which  is  qualitatively  the  same  as  the 
stress-strain  relation. 

Before  solving  a  complicated  plane  strain  problem,  we  next 
check  the  algorithm's  ability  to  simulate  rotational  rigid-body 
motions.  This  is  done  by  fixing  a  corner  of  the  stressed  element 
and  prescribing  the  essential  boundary  conditions  at  each  incremental 
step  as  in  TAYLOR  and  BECKER  [1983].  Suppose  that  a  block  element, 
such  as  that  in  Fig.  5.4,  is  subjected  to  an  elasto-plastic  defor¬ 
mation  according  to  the  following  program:  We  prescribe  0.5  percent 
engineering  strain  with  a  0.5  x  10  per  second  rate.  Next  the 
block  is  rotated  with  prescribed  incremental  rotation  angles  while 
maintaining  the  preloaded  deformation.  The  results  are  listed  in 
Table  5.1  with  increments  9°,  10°,  15°,  30°,  and  45°  to  make  90° 
rotation.  If  the  computed  stress  rates  are  appropriately  objective, 
the  stress  should  not  change  during  these  rotations.  We  observe 
errors  in  the  stress  of  only  about  0.3  percent  and  note  that  there 
is  little  error  between  a  large  step  size  (45°)  and  a  small  (9°), 
as  expected  from  the  way  we  defined  deformation  measures  and 
rotational ly-invariant  stresses. 


strain 


S,  =  1  6  1.3  N/  mm 


_ J  J2  =  80.7  N/  mm 


(a) 


(b) 


Calculated  rigid  body  motion  with  30°  increments, 
a)  loading  before  rotation 

progressive  configurations 


b) 


The  final  example  is  a  head  forming  problem  in  plane  strain. 

A  4  by  5  unit  rectangular  billet,  which  is  confined  at  the  lower 
boundary  is  loaded  at  the  upper  part  without  friction.  Incremental 
displacements  are  prescribed  at  the  five  nodes  in  the  top  of  the 
billet.  Computed,  progressive  deformed  shapes,  and  stress 
contours  are  shown  inFigs.  (5.6)  and  (5.7).  Figs.  (5.5)  and  (5.8) 
show  the  undeformed  and  the  deformed  Lagrangian  finite  element  mesh 
The  residual  stress  contour  is  shown  in  the  left  part  of  the 
deformed  configuration. 

Throughout  the  finite  element  computations,  the  convergence 
at  each  incremental  step  was  checked  by  calculating  the  maximum 
relative  error  of  successive  incremented  displacements.  The 
relative  error  is  computed  as  the  ratio  of  the  correction  between 
iterations  to  the  first  solutions  (incremental  displacements)  of 
the  incremental  step.  A  range  of  tolerances  was  set  as  0.01-1 
percent,  depending  on  step  size.  Generally,  the  convergence  was 
achieved  in  two  iterations  except  where  severe  changes  in  the 
deformation  from  elastic  to  plastic  states  are  experienced. 


CHAPTER  VI 


CONCLUDING  REMARKS 

6.1  General 

The  classical  theory  of  plasticity  evolved  from  simple  obser¬ 
vations  that  irrecoverable  deformations  of  many  materials  may 
result  from  a  simple  load  cycle.  Characteristic  of  the  classical 
theory  and  the  developments  which  arose  from  it  over  the  last 
half-century  are  notions  of  yield  surfaces,  hardening  laws,  and 
a  rather  simple  kinematics  in  which  the  total  strain  is  the  sum 
of  so-called  elastic  and  plastic  parts.  Frequently  in  modern 
engineering  and  manufacturing,  situations  are  encountered  in  which 
a  more  sophisticated  constitutive  theory  and  kinematics  is  needed 
to  describe  elastoplastic  behavior.  This  need,  and  also  the 
natural  evolution  of  plasticity  theory  as  a  part  of  mechanics, 
has  led  to  some  fundamental  questions  of  the  mechanical  and  math¬ 
ematical  foundations  of  plasticity.  For  example,  how  does  one 
develop  a  theory  of  elastoplasticity  within  the  framework  of 
modern  continuum  thermomechanics?  How  does  one  obtain  a  set  of 
constitutive  equations  which  are  valid  for  wide  ranges  of  defor¬ 
mations?  What  are  the  deficiencies  of  existing  kinematical  measures 
and  is  it  possible  to  devise  new  kinematical  measures  appropriate 
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for  finite  elastoplastic  deformation?  Can  such  general  theory 
be  put  to  use  to  model  practical  problems  and  is  it  possible  to 
devise  numerical  schemes  for  implementation  of  such  a  broad  theory 
In  this  report,  an  attempt  is  made  to  resolve  some  of  these 
questions.  Some  success  in  each  area  of  inquiry  has  been 
attained. 

6.2  Summary  of  Results 

According  to  the  general  aim  of  this  report,  the  following 
specific  results  were  obtained. 

1.  Materials  of  Type  N 

A  new  and  general  theory  of  finite  elastoplasticity  has  been 
developed  which  has  the  following  features: 

a.  Involves  two  potential  functionals,  one  being  the 
free  energy  functional  <p  and  the  other  the  general 
flow  potential  < l>. 

b.  Does  not  necessarily  require  the  specification  of 
a  yield  surface. 

c.  Does  not  require  that  ip  be  either  convex  or  differ¬ 
entiable  . 

d.  Involves  an  internal  state  variable  a. 


e.  Reduces  to  classical  theories  as  special  cases. 


f.  Is  not  exclusively  a  rate-type  theory,  as  the  history 
of  microstructural  changes  is  governed  by  an  evolution 
equation  for  a. 

The  material  characterized  by  these  features  is  referred  to 
as  a  "materials  of  type  N"  since  the  framework  generalizes  the 
notion  of  normality  in  plasticity. 

2.  Kinematics  of  Finite  Elastoplasticity 

A  new  kinematical  description  of  notion,  valid  for  finite 
elastic  and  plastic  deformation  and  which  differs  from  the  theories 
of  Lee  and  Nemat-Nasser  has  been  derived.  Its  features  are  as 
follows: 

a.  The  effects  of  rotation  are  separated  from  pure 
stretching. 

b.  Formalizes  the  fact  that  the  position  vector  (or 
displacement  vector)  in  the  intermediate  configuration 
may  not  be  continuous. 

c.  Results  in  a  correct  decomposition  of  deformation 
rate  into  an  "elastic"  and  a  "plastic"  part. 

3.  Continuum  Thermomechanics  and  Analysis 

a.  Analysis  via  Coleman-Nol 1  thermodynamic  methods  are 
used  to  deduce  restrictions  on  the  constitutive 


equations . 


b. 


It  has  been  shown  that  the  theory  results  in  a  fully- 
determined  system  of  equations  governing  the  thermo¬ 
mechanical  behavior  of  materials  of  type  N. 

c.  Isotropic  representations  of  flow  potential  and  the 
free  energy  for  the  elastoplastic  material  of 
materials  of  type  N  has  been  developed  and  studied. 

4.  Examples  of  Materials  of  Type  N 

a.  Classical  plasticity  theories  are  shown  to  be  special 
simple  cases  of  a  materials  of  type  N. 

b.  An  analyses  of  some  newly  developed  constitutive 
equations,  which  do  not  involve  yield  functions,  is 
given. 

c.  Existing  constitutive  equations  for  various  plasticity 
theories  can  be  shown  to  fit  within  the  theory  of 
materials  of  type  N  upon  appropriate  modification. 

5.  Finite  Element  Computations 

a.  Chooses  a  modified  version  of  the  Bodner  and  Partom 
equations  for  finite  element  approximation. 

b.  Numerical  algorithms  and  a  code  are  developed  for 
solving  a  total  Lagrangian  formulation  of  the  general 
theory . 

c.  A  test  problem  of  rotation  of  a  prestressed  block 
and  head  forming  problem  are  solved. 


6.3  Future  Research 


1.  Further  study  is  needed  to  verify  that  other  useful 
constitutive  equations,  e.g.,  by  HART  [1970,  1976],  by  MILLER 
[1976]  and  anisotropic  version  of  BODNER  [1985],  can  be  recast 
in  a  form  covered  by  the  theory  of  materials  of  type  N.  It  is 
anticipated  that  such  modifications  are  possible  and  their  con¬ 
struction  would  be  very  helpful. 

2.  To  predict  the  useful  lifetime  or  realibility  of  perfor¬ 
mance  of  engineering  material,  failure  or  damage  accumulation  laws 
need  to  be  developed.  It  is  expected  that  it  may  be  possible  to 
incorporate  these  effects  in  a  theory  through  the  introduction  of 
appropriate  internal  state  variables.  One  possible  example  of 
this  may  be  the  work  by  LEMAITRE  [1984]. 

3.  Further  development  of  a  modular  computer  code  is  needed 
which  can  simulate  more  realistic  problems  and  may  include  temper¬ 
ature  effects  and  appropriate  frictional  contact  conditions, 
possibly,  in  a  three-dimensional  setting. 
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APPENDIX 

ELEMENTS  OF  CONVEX  ANALYSIS 

We  shall  provide  here  a  brief  summary  of  some  of  the  concepts 
of  convex  optimization  theory  which  are  prerequisite  to  the  ideas 
described  in  the  body  of  the  report.  For  more  detailed 
accounts,  the  books  of  EKELAND  and  TEMAN  [1976]  or  ROCKAFELLAR 
(1970,  1979]  or  the  recent  text  of  ODEN  [1985]  can  be  consulted. 

We  begin  by  introducing  the  following  notations: 

HI  =  the  extended  real  numbers;  if  1R  is  the  real  number 
system,  IR  =  IRU{±°°} 

U,V  =  topological  vector  spaces 

U*,V*  =  topological  dual  spaces  of  U  and  V  respectively 
(.  ..)v»  <•»•> jj  =  duality  pairing  on  V*  x  V  and  U*  x  U, 
respectively;  i.e.  if  v*  f  V*  and  v  £  V,  then  v*(v)  h 
<v*,v>v,  etc. 

It  is  worthwhile  to  recall  the  definition  of  the  limit-superior 
(lim  sup)  and  the  limit-inferior  (lim  inf)  of  sequences  of  real 
numbers,  extended  real-valued  functions  of  sequences,  and  sequences 
of  sets  in  a  topological  vector  space  V. 

tin  6up/tuv  ini. 

'  For  {a^l  a  sequence  of  real  numbers 


lim  sup  a  =  inf  sup  a 
n  n 

n  ♦  «  N  €  IN  n>N 


lim  inf  a  *  sup  inf  a 
n  n 


N  €  IN  n>N 


For  f:  1R  *  1R 

lim  sup  f(x')  =  inf  sup  f(x') 

x'  -  x  6>  0  0  <  |x'-x  |  <6 


lim  inf  f(x')  =  sup  inf  f(x') 

x'  ♦  x  6>  0  0  <  |x'-x|  <6 


For  {A^  }  a  sequence  of  subsets  of  the  underlying  set  of 
topological  space  V, 


lim  sup  An  -  (\J\) 


m“ 1  n=m 


lim  inf  A n=\J  Am) 


n  -►  oo 


m=l  n*m 


For  example,  if  {x^  }  is  a  sequence  of  real  numbers  which  converges 
to  x  and  f  :  1R,  LIM  SUP  f(xn)  is  the  supremum  of  all  cluster 
points  of  f  at  x,  as  indicated  in  Fig.  A.l  (with  an  analogous 
interpretation  for  lim  inf). 

The  concept  can  also  be  applied  to  multifunctions  from  one 
topological  vector  space  to  another.  Indeed,  if  T  :  U  ♦  V  (with 


r(u)  a  subset  of  V  for  each  vector  u  t  U),  then 


lim  sup  r(u’  )  «  i  \  [\J  (r(u')  ♦  A)1 

u’-u  A  6N(0)  BfW(u)  u'fB 


and 


lim  sup  T(u' )  = 
u '  ♦  u 


A 


(A. 5) 


A6W(0)  BfW(u)  u'  £  B 


(r(u’)  ♦  A): 


where  b!(0)  and  N(u)  are  collections  of  neighborhoods  of  0,u, 
respectively. 

Lim  iap  inf/ Lim  inf  6up 

In  addition  to  the  notion  of  limit  superior  and  limit  inferior, 
it  is  convenient  to  introduce  the  concepts  of  lim  sup  inf  and  lim 
inf  sup  introduced  by  ROCKAFELLAR  [1980]. 

Let  F  be  an  extended  real-valued  function  from  U  x  V  into  IR, 

let  u'  *  u  in  U  and  v' ♦  v  in  V.  Then  we  define 

A  (A‘6) 

Lim  sup  inf  F(u'v')  *  sup  inf  sup  inf  F(u’,v') 

u'  ♦  u,v'  ♦  v  B  6  W(v)  A  f  W(u)  u*  A  v'(B 

Likewise , 

A  (A'7) 

Lim  inf  sup  F(u',v’)  *  inf  sup  inf  sup  F(u',v') 

u'  *  u ,  v '  ♦  v  (v)aN(u)u'€Av'€b 

Similarly,  Limsupsup  and  Lim  inf  inf  can  be  defined  in  an  analogous 
way . 

The  meaning  of  these  operations  can  be  more  easily  understood 

2 

in  the  case  of  a  real-valued  function  F  defined  on  IR  ,  such  as 


the  discontinuity  at  the  origin  shown  in  Fig.  A. 2.  To  compute 


Lim  inf  inf  F(x',y),  for  example,  we  compute  lim  inf  F(x',y)  for 
x'  *  0,  y'  *  0  x'  -►  0 


a  fixed  y'.  This  gives  the  a  function  of  y*  which  has  as  its  graph 

the  curve  AC  U  EH.  The  lim  inf  of  this  function  is  the  point  E, 

y '  ♦  0 

denoted  F  in  the  figure.  Similarly  lim  sup  F(x',y')  for  fixed  y', 
1  x*  -  0 

is  the  curve,  HD  U  BA  and  lim  inf  of  this  curve  is  the  point  D, 

y'  +  0 

denoted  F^  in  the  figure.  In  summary  for  this  example, 

Lim  inf  inf  F(x',y')  *  F. 
x 1  0  y '  ♦  0 

Lim  sup  inf  F(x',y’>  =  F_ 
x'  -  0  y'  -  0 

Lim  inf  sup  F(x',y')  *  F. 
x'  -  0  y'  -  0 

Lim  sup  sup  F(x',y')  *  F, 
x'-  0  y'  -  0 


J 


V 


"W  " 


-  1 


F7“  V  • 
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